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ABSTRACT 


ee GollGATION OF THE STABILITY OF A SYSTEM OF 
TWO SHIPS EMPLOYING AUTOMATIC CONTROL WHILE ON 
PARALLEL COURSES IN CLOSE PROXIMITY 


Charles Natale Calvano 


submitted to the Department of Naval Architecture and Marine 
Engimeering on 21 May 1970 in partial fulfillment of the 
requirements for the degrees of Master of Science in Ocean 
Engineering and Naval Engineer. 


An introduction outlining the development of ship 
Operations in close proximity for underway replenishment is 
followed by a survey of the literature Ceding sy seme the 
measurement of interaction effects between ships under such 
circumstances. A brief outline of the theoretical treatment 
of ship motions in open water is given and followed by the 
modifications to the theory required for treatment of the two 
ship system. The applicable equations are written for calm 
Water and suitably linearized. The techni ques (Of meneame unc 
the differential equations as algebraic equal@ons, amemot 
solving the equations by determinant evaluation are described. 


Extension of available interaction and open water data 
FOr two different ships to a full data set for a single, 
Meometed Ship 25 carried out, and the numerical values of 
resulting hydrodynamic derivatives are given. The numerical 
SsOlution of the equations and the extraction of Saeco dale 
roots are then described and results for 9] typical trial 
data sets are listed. A discussion of the results obtained 
is followed by an evoluation of the correctness of the 
Solution method employed. 


The computer programs used cthrougheut aremde-cr mea an 


detail in four appendices. 
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Chapter J. 


INTRODUCTION 


ae Packground-Replenishment at Sea 

In their efforts to achieve a degree of independence from 
Hand bases, the navies of the world have, for decades, prac- 
ticed various methods of underway refueling and replenishment 
of their ships. As techniques were imoroved and experience 
was gained, these operations progressed from the status of 
mreecquent, makeshift affairs, to that of routine evolutions. 
Giese replenishment methods have ranged from actually lashing 
together the supply ship and receiving ship while passing 
coal and provisions from one to the other, (one assumes this 
was attemoted in only the very calmest of seas, or in shel- 
mered roOadsteads), to the transfer of cargo via "vertical 
meorenvshment” (WERTREP) helicopters. Other methods ant ome 
over the years have included the astern refueling method and 
the alongside method for transferring fuel and cargo. 

In the former, the receiving ship recovers a fuel hose 
trailed in the water by the tanker Vara takes on fuel while 
Steaming in the tanker's wake. In the alongside method, the 
receiving ship overtakes the supply ship on a parallel ccurse, 
laterally separated by froin 40 to 180 feet. The transfer is 
effected while the ships remain abreast at this distance. 

The alongside method is the standard of the U.S. Navy, 


(ceuqmem@ceds by VERTREP for cargo transfer, though not, of 


course, for refueling), and such replenishment onerations 





meee oecome almost a way of life for U.S. combatant ships. 
Beress, and Until, all U.S. combatant ships become nuclear 
meeted, this method will probably ere in wide use. 

Of course there are other situations than replenishment 
when two ships might wish to operate in close alongside prox- 
imity while underway. To generalize the treatment herein, 
mee forces and moments arising from the tension in the 
replenishment rig lines is not included. Thus the situation 
occurring whenever two ships must steam alongside each other 


GImemr@sce Quarters is that which is treated in this thesis. 


mi  eome Problems 

Although alongside replenishment has become commonplace, 
this type of close-in steaming has not ceased to be a dangerous 
undertaking. There are, of course, dangers arising from ea 
ack of freedom to maneuver to avoid other shipping, or, in 
tame of war, the actions of an enemy. But the ever-present 
source of danger to two ships steaming close together on 
Perallel courses is this very proximity of the ships, causing 
Seelweremeenstitute a collision hazard for the other. 

While alongside, the ships must make almost constant 
rudder and cngine speed adjustments to maintain position. 
Because the ships are close, the presence of each affects the 
forces and moments on the other through the influence of each 
on the pressure field and flow patterns around the other. 
(See Figure I-1). Thus, neither ship will respond as it does 


MMC vwetler where these influences are not présemt. And 
4 a 
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Peessture Field about a Shap in Open Water 
(From Reference A-1) 


FIGURE I-1 


as the ships alter their relative positions, the direction and 


magnitude of these influences are continuously altered. 


The repeated corrections required over a period of up to 


several hours, introduce hundreds of opportunities for human 


Gmeor Or momentary lack of attention to give rise to damaging 


conseauences in the form of a collision. Any method or device 


which could serve to reduce the probability of such errors 


Gecurtine would serve to lessen a primary hazard for two ships 


Gretatane in this way. 


va 





meee SULOMatic Control] 

The antroduction of automatic rudder angle and engine 
Speed control to the alongside operation is an attractive 
prospect. Such a control system would employ sensors to 
measure longitudinal and lateral separation of reference 
points on the two ships, (and their rates of change), and 
compare these values to desired equilibrium values. It would 
then generate rudder and engine orders to reduce the deviations 
to acceptable Peonitudes. such a system could be a boon to 
conning officers and helmsmen by relieving them of much of the 
tedium of maintaining position alongside. 

A theoretical investigation of such a system is the 
Subject of this thesis. The problem will be approached from 
a mathematical standpoint as a two-body system of ship motions. 
The hardware required to make it operable will not be examined 
or described. It will be assumed that a sensor, (doppler sonar, 
for instance), 1S available to continually measure lateral and 
longitudinal Ship separations and their rates of change. 

Once the deviation from desired ship separation and heading 
1s known, an automatic control system can order rudder and 
engine speed changes to reduce the deviations. The various 
constants associated with a theoretical control system meeting 
this description will be a part of the mathematical model of 
Etemiwo Ship system. The stability of the system wili be 
investigated, to determine if it appears feasible to employ 
Such an automatic control system to keep the ships in position 


alongside. (By stability is meant the tendency of the two 


V2 





Mies yotem tO return to an equilibrium condition after 

receiving a small disturbance from that equilibrium condition. 
[Meecthe System ‘returns, or tends to return, to the equilibrium 
Seereronr when the disturbance is removed, then it is sta@ble”. 


m3 1*) 


Mmeeiastorical Aspects--The Literature 

Peaeelayilor. 

As far as the author has been able to determine, the first 
Systematic attempt to measure interaction phenomena for two 
Ships on:parallel courses in open water was made by David W. 
Taylor, the ubiquitous Naval Constructor. He conducted a 
series of tests utilizing four models, each of 3,000 pounds 
displacement, at the Model Basin in Washington, aia) 1900. bee) 

Taylox's tests were conducted by towing the models in 
pairs, either east of one another, oO With ene meee 
fixed distance ahead of the other. The models were not self- 
propelled, and did not have propellers affixed. They were 
towed under towing carriages, with their bows and sterns 
prevented from translating laterally, but free to move vertic- 
aiy, Changing draft and trim while being towed. Water depths 
were “many times the draught of the models". Attraction and 
repulsion Pome were measured at the attachment points near 
tmembow and stern of each model. Table I-l1 Jists the eharac— 
teristics of the models employed and Figure I-2 is a typical 
plot of the results Taylor obtained. 


*fetter-number citations in brackets refer to the references 
listed at the end of the body of this thesis. 
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Preeenauccting his tests, Taylor found that he hed 


Ger cculty towing the models on a straight course. In view of 


the fact that towing carriages were used, this is hard to 


understand, but was explained no further. As a result of this 


adifficulty, he felt that the forces measured could not be 


regarded as highly accurate, but felt that they showed the 


tendencies and general nature of the interaction, or, to use 


his term, 


"suction", phenomena. 


The sequence of events, when vessel B overtakes vessel A 


on a parallel course, seemed to be about as follows: 


was 


als 


As B's bow overlaps A's stern, a small repulsion 
force is felt at both B's how and stern. 

AS B continues to overtake A, B's bow is attracted, 
while its stern is repelled, causing B to tend to 
iam ENGOs A. 

As B draws abreast of A, both B's bow and stern 

are drawn in toward A. 

As B pulls up farther, the force on its bow switches 


to repulsion, while its stern is still attracted 


to A. 


Peo bp draws farther forward, to where ts sterm ae 
approximately abreast of A's midlength, both B's 


Maw dame stern are repelled: 


Tayler noted further, that the occurrence of titese «va. 


somewhat more consistent with the "finer" models, while 


tremeesnlts with the “full" mdédels tended to be”™erratic. The 


Cmeecte of 


the “suction" forces on the models' resistance was 





MmoteMeaeerca, but Taylor did note that the side forces "varied 
With Speed as the resistance of the model"; that is the inter- 
action forces increased as the square of the velocity. 

It is evident from the discussions accompanying Taylor's 
paper that there was disagreement within the Naval Architect 
community concerning the significance and, in fact, even the 


existence, of interaction forces for full scale ships. 


B. Gibson and Thompson. 

In 1913, Gibson and Thompson conducted a series of full 
scale tests with small steam launches in the Thames River 
near London, to verify the existence, and measure the effect 
of interaction phenomena. [B-2] They did note effects 
Samilar in type to those described by Taylor, but due to 
instrumentation difficulties, their report contains no 
Quantitative information of a useful nature in considering 


iacberaction effects for ships. 


feomeast ous Collision between the liner QUBEN Tiehveand 
the destroyer H.M.S. CURACOA in 1942 led to a series of 
experiments reported by A.M. Robb in 1949. [A-8] In these 
tests, the models used were of these two ships and made to a 
1356 scale. This resulted in the QUEEN MARY model being 
“rather more than 18 feet" in length, and that of the CURACOA 
"hatmer more than 8 feet". In these experiments, both models 


were se]. f~propelled, with first one, then the other, being 


Ceomemtratned tO run On a straight course while the other passea 
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On a nearly paraliel course. Robb again confirmed the exist- 
ence of significant interaction effects drawing the smaller 
eaip into the larger, in some circumstances, and his article 
1s accompanied by several sketches showing the trajectories of 
the ships for various speeds and lateral separations. The 
Peeults presented are largely qualitative in nature, and do 
not disagree with those obtained by Taylor 40 years earlier. 
But again, in the discussions of Robb's paper, the fact 
that the very existence of interaction effects was doubted 
in some quarters was brought out. [B-5]) Professor Prohaska, 
in his discussion, stated, "The question of interaction between 
passing vessels has always given rise to interesting exchanges 
Grmmoeriion..... a great many experts, amongst whom are numerous 
seafaring peovle with great experience, deny the existence of 
interaction, or at least think it without importance in 
unrestricted water." (Prohaska, however, did not share this 


View.) 


BP eoericted Water Studies. 

Over the years, there have been ee other experiments 
and theoretical developments concerning the interaction between 
a Ship and a stationary, constant-geometry object, such as a 
Gdmerakin Or fbenk. References A-6, A~LO, B-6,,Bsi0,_ B-11, 
Pee lo ana B-18 all deal with this phenomenon. A ship, 
Demiena large object, can appear, in its effects, as a wail 
would appear to another ship passing near. But of course, 
thewowete ry Of the ship-object is different than when che object 


1s a wall. And, also, neither ship is immovable, as is the wall. 
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Most of the investigations concerning interaction in 
@enalseor channels include shallow water conditions. The 
shallow water tends to emvohasize the interaction effects caused 
by thewwall. This investigation, though, 1s concerned with the 


interaction between ships in deep water. 


E. Newton. 

The single study of interaction effects which prcved to 
be of greatest use to the author in attempting to obtain 
quantitative information concerning the interaction forces and 
moments was that listed as reference A-l, by R.N. Newton and 
reported at the First Symposium on Ship Maneuverability 
Cemaucted at the David Taylor Model Basin in 1960. 

Newton reported on experiments conducted at the Admiralty 
Exveriment Works (AEW) in England, between 1946 and 1948 and 
which were inspired by the growth of alongside replenishment 
experienced uleene tne war. In his paper, Newton acknowledges 
the pioneering work of Taylor in this field and laments the 
near~-total lack of work in this area in the years between 
Taylor's tests and the AEW efforts, although he noted the 
wealth of information available to him concerning interaction 
effects in restricted channels and shallow water. 

Newton's experiments were with models of a battleshin, 
H.M.S. KING GEORGE V and a revlenishment supply shin, R.F.A. 
OLNA. (See Table I-2 for ship characteristics.) The tendencies 
exhibited by the models were confirmed by full scale trials at 


SG. 
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a ee ets + ae Oo eee = ep ee er ee 


OLNA 
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Pengen on W.L., £t. . 547 

Beam, ft. 70 

Diem: , ft. a0 

Displacement, tons 25757 0 

Block Coefficient i/14 
Corresponding water depth 75 fathoms 
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Tests were conducted with both towed and self-propelled 
versions of the models. The models were allowed complete 
vertical freedom, but constrained from Jateral motion. Forces 
and moments were measured at longitudinal sevarations of the 
[Mee sections Of the ships from 600 ©t. aft of the abmeast 
Poetenen to 600 £t. forward of this position, at two lateral 
Separations; 50 oie UOOM Eta Sade to side aqistance. Wicnhmehe 
Ships in-the abreast position, the side to side distance was 
varied from 25 to 140 ft. (All distances are the corresponding 
full scale values.) 

ieeiteesce t= and L-4 showy the longitudinal variation of 
the forces and moments experienced by OLNA at the two lateral 
separations and Figure I-5 shows the lateral variation with 
the ships abreast. These data formed the basis for quantitative 
interaction forces and moments used in this investigation, as 
Oi wes NOt radically different in hull characteristics from 
the MARINER class vessel used for this study. In the figubles. 


OLNA is the ship being overtaken. The agreement with the 


~ 


Paoweeeos Oocalned by Taylor is apparent. 
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Hee econ lesgsenvand Chislett. 

In making a theoretical investigation such as this one, 
it is necessary to have realistic tines mor thebhivdrodynamic 
derivatives employed in the mathematical model. (The use of 
these derivatives, and the mathematical model are discussed in 
Cieweter Ll, THEORY.) For two ships steaming in proximity, 
values of both the open water and interaction derivatives are 
needed, as forces and moments which are always present for a 
ship in open water will be present along with those caused by 
the interaction phenomena. 

Thus, in this investigation it was necessary to find 
SOmecce ten aaypaecal values of both kinds of hydrodynamic 
derivatives. As mentioned earlier, Newton's work with KING 
GEORGE V and OLNA was the only available source of quantitative 
Mimecact ten innonmationy It was, therefore, necessary voce 
values of the open water derivatives for a ship similar to 
either KING GEORGE V or OLNA. Strom-Tejsen and Chislett in 
reference A-5 determined most of the open water hydrodynamic 
derivatives of interest for a MARINER class merchant ship. 

(eee, taole 1-3 for MARINER characteristics.) In reference A--4, 
strom-Tejsen used these derivatives to theoretically predict 
the open water maneuvering behavior of a MARINER, and 
subsequent sea trials proved these values to be quite accurate. 
Fortunately, OLNA and the MARINER class hulls are sufficiently 
Similar to expect the hydrodynamic derivatives of one to 
reasonably resemble those of the other. Therefore, in the 


development of the results of this thesis, the open water 
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derivatives of the MARINER were used, and, to account for 
interaction effects, derivatives derived from the interaction 
force and moment data of the OLNA experiment were added to 
these to obtain the entire list of needed derivatives. MThus, 
the ship type investigated is a somewhat modified MARINER 


class ship. 





MARINER CWecacCucr Hse les 





bemgth, “rt. 





Beam, ft. 






Diwert, rt. 





Displacement, tons 






Pilock CoelLfEicient Cy 


Piacsmatic Coefet. Cy 






Midships Section Coeff. oe 





me 


TABLE 1-3 
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Ciapecr, If 


THEORY 


a ee eee 


Bo Background 

The derivation of a mathematical model representing the 
steering and maneuvering of a surface ship in open water is 
given by Abkowitz in reference A-3. In its general form, the 
model accounts for non-linear, as well as linear, effects, 
and in his computerized prediction of open water, surface ship 
Maneuvers, Strom-Tejsen uses the Abkowitz model with up to 
Smmcoraer terms. [A-4] 

This mathenatical model, in linearized form, serves as 
the basis for the analysis made in this thesis, with appro-~ 
priate extension of the model to include the interaction 
effects which are not present in open water, but which are 
GEmecrrane FOr two Ships on parallel courses in close proximity. 
The model used herein also is modified to represent the two- 
ship system, rather than the motions of a single ship, and 
accounts for the fact that both ships are free to move in the 
horizontal plane. (For a more complete description than the 
outline which follows, of the derivation of the open water 


model, reference A-3 is suggested). 


II. Reference Axes 

Pew beletrence axes used constitute an orthogonal, right- 
handed system which takes advantage of the existence of a 
vertical plane of symmetry through the centerline of the ship. 


} 
Peet octiOn Of the axis system illustrated introduces 
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Simplifications in the equations of motion by virtue of the 
fact that the plane of symmetry of each ship is parallel to 


the principal axes of inertia of the ship. 


S TED. 


oa 








vty : 


(See Table II-1 for definitions of all symbols used 


m1 this chapter) 


Reference Axis System 


PreuRe -ii—} 


Pieter ons Of Motion for a Ship Moving in the Horizontal 
Plane 
A general form of the equations of motion for a body 
which is ProtmecieS move in all the six degrees of freedom 
iS given by Abkowitz. 
Normally, when dealing with steering and maneuvering of 
surface ships in open water, the primary motions are considered 


€06 take place in the horizontal plane, and vertical motions are 
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neglected. As this is also the case for the two-body system 
Of interest, where only horizontal motions are of concern, 
meme equation in Z, the vertical force equation, can be 

Soo oead. FUrther, choice of the axis system parallel to the 
plane of symmetry of the ships eliminates the cross terms of 
M@erecia, and, if 1t iS assumed that the center of gravity of 
each ship is in its centerline plane, €@s, indeed, can be 
expected), the quantity Ve is zero for both ships. The 

usual assumption that rolling and pitching motions, if not 
extreme, have little influence on steering and maneuvering is 
also made , Presi the equations for K and M ae be neglected, 
and making the variables p and Cee Oy. (See Table II-1 for 
definitions). Thewequations of motzroen forea singlessurface 


shio moving in the horizontal plane thus become: 


a = m[Q-rv-x 5°] 

ba Lo a 

= ie Ti+mx, (v+ru) (2=)) 
(where the notation v = dv/dt, is used) 


iiewlc teaiadnemsctacs Of these equations represent the sum 
Cue LoOuces an tie xX and y directions and all moments about 
Eiem~enaitas, respectively, .The raght~hand sides are the 


dynamic response terms associated with the ship and its motions. 


IV. Taylor Expansion of Forces and Moments 
ihe LOrces#ama momentseon the hefét-hand sides»of the 
equations of motion can be expressed as functions of the 


properties of the body, properties of the fluid and properties 
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TABLE TI-1 


Definitions of Symbols 
Reference coordinates fixed with respect to 


earth. 
Reference coordinates fixed in ship. 


Melloctieves an Glrection Of x, y and 2 axes, 


respectively. 


Roll, pitch and yaw angles about x, y and z 


axes, respectively. 


Angular velocities about x, y and Z axes, 


respectively. (p=¢, g=6, r=). 


Rudder deflection in radians (port rudder is 


positive). 
Propeller speed in revolutions per second. 


SUNG One sOsces SacuIng On Ship in x, vy Zz 


directions, respectively. 


Silo OLeMOMentsS acting On Ship about x, y, 2 


axes, respectively. 


Longitudinal separation of ships measured 
from Ship A midships section to Ship B 


midships section. 


Lateral separation of ships measured from 


Ship A centerline to Ship B centerline. 


Typical of notation used for hydrodynamic 
derivatives. In this case: d3Y/du for Ship B, 
measured when all variables but u are at 


equilibrium values. 


Control system time lags associated with rudder 


and propeller response, respectively. (Subscripts 


). 
Le 


A Or B indicate applicable ship, as At 
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COmerOlMewoeem Constants, 1=1 to 6, for Ship aA 
or B. Constants associated with control system 


response to deviations from equilibrium, as 


§. =k, Ay. 
A oe 


Distanee Ol Center Of Gravity of ship, (G), 
Pac iewerers OL1G lye (0), in x, y, z2 directions 


EGS Weertvely. 
Mass of Ship A or B in units of (lb sec’/£t) 


Hite OweisatOr d/ dts, asec a= ou/ ot 
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of the motion. When considering specific hull forms, and 
using generally accepted scaling laws, the forces and moments 
may be considered functions of cnly the motion and orientation 
parameters of the bodies. (When dealing with automatic 
control of rudder angle and propeller speed, they are also 
functions of rudder deflection, 6, and propeller speed, n.) 

The reference system for the two-ship system at 
equilibrium and an illustration of the variables involved 
are shown in Figure II-2. 

When considering motion in an unrestricted horizontal 
Plater it is Clear that no Eorces Or moments arise from a 
change in ship location with respect to a fixed geographical 
Mr ence point, allowing Xo and we to be discarded as. 
variables of interest. Similarly, the eguilibrium yaw angles, 
CCWad5 and Wp),), ox drift angles, can be drepped from the 
list of independent variables because forces and moments 
eisai L£rom a drift angle are accounted for by the terms in 
uand v. (Tf there were no drift angle, u would coincide with 
velocity made good in the desired direction, and v would be 
zero. Thus the forces and moments that Wary and ad, cause 
to be introduced are taken up in the terms in wu andv. The 
inter-relation of drift angle and the independent variables, 
Weanedev, can be seen frem tan(drift-angle) = v/u). 

It is.also possible to drop a, O, 8 and 6 as independent 
variables because of the simple eeteonencerie relations 
between them and v, Vv, u and &. The usual assumption that 


forces and moments produced by ae § and n are negligible, in 
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fermenor sip Control, 1s also made. (These forces and 
moments are not necessarily negligible when considering forces 
On the rudder stock, propeller blades, etc. But they do have 
a negligible effect on shiv motion. (See reference A-4, 

pg- 7).) 

As Abkowitz has shown, the left-hand sides of the 
equations of motion, comprising the functions describing the 
forces and moments acting on the ship, can be developed into 
Wweerul form for ay sis iby Wse "Gr Ene Way lor series vex eansien 
for a function of several variables. The force or moment 


¢ eet 


caused by a change in any variable is expressed as the product 
——_—eeearr en a ae ~ - 





of the derivative of the force with respect to that variable, 
(with all other variables at equilibrium values), and the 
ane in the variable. In the linearized form of the 
equations, only the rear terms are retained, though the 
expressions can be written to any desired degree of accuracy, 
by retaining sufficient terms in the expansion. (For example, 
oe 4 force aenens by a change in u is written as (dy/du)x (Au) 
where (dY/du) is the value of (dY/du) with all other variables 
at their equilibrium values. The derivatives, such as (dY/du) 
are Bi Pig EARN CeLiIvavivics ana. Can OC Cxperrmnen tansy 
Getemnamed: ) 

Employing standard notation, where ‘an represents (dY/du) 
at equilibrium, the Y equations, for example, for the two 
ships, retaining only the linear terms of the Taylor expansion, 


become: 
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Ship A: 


a) Y, Aa +¥ AB+Y,., Aud. OU,tY Aviat... . 


A Bn / A A A 
a” 
-+Y « AS,tY, Any == el ry] 


A A A 


Een eno 
O 


“- + er 
: AB pe Aunty. Au, e. Avit+ 


Pp B B a 


Ty Ye NOY) 6A C= [vi+r_u_+x fr] (2-2) 


eewton es | ob BBB g. B 


The left-hand sides of the X and N equations are similarly 


_expanded for both ships. 


(eee atizing the Dynamic Response Terms 

Since only linear terms have been retained in the 
expanded force and moment expressions on the left-hand sides 
of -the equations, a Similar linearization of the dynamic 
response terms on the right-hand sides must be made. 


Poreall variables, the following holds true: 


Variable = (variable) | + Avariable (2-3) 


nN t! 


where the subscript "o" refers to the equilibrium value. 
for all time derivatives of the variables, (v, v, ere. 
the equilibrium value is zero, as the equilibrium condition 


is unchanging in time. Thus for these: 


Tar iabLles= (variable) | +ixencaeeoukte = Awarvable (2-4 ) 


I 
a 
i 
G 


(ec. U = ur) 


Makang use of (2-3) and (2-4), and by neglecting second 


ao 





mage eee ris CONsascting Of products of small quantities, such 
as (Ar)(Au), the linearized right-hand sides of the Y equations 
become: 
Y=m{[v + uor + ee 3 (2-5) 
Similar developments are applicable to the X and N 


equations and are carried out in detail in [A-3]. 


VI. Linearized Equations of Motion 
PoamiagelitrO accoun. all developments Outlined thus far, 


the linearized equations of motion become: 


Ship A: 
a) X AatX, AB+X. Au_+X. &.+XK Av.+X. V 
On Ba | Uy A Uy A Vn ox Vn A 
+X* WitX. botX, AS +X. An. = mu 
va A uN A 6 5 a hy A A~A 


—-. 
b) a AatY AB+ x Au, +Y.. uty Av,tY. V 


A ee, A A A A 


Yh vat oe Oe, — Me to 


Va A A A 


c) N. AotN, AB+N. Au_4+N. U_L+N Av_tN. V 
o u A -u V 


A Pa A pe Wigan Ya 2 
+N° tN. §.4N, AS +N. An. = I jotm x (v.4u_ p,) 
eee tA Sos Aon, A z,/A" A°g, A '0,'A 
dq) 
e) Comracesmonarniag law@ations cor Ship B 
i.) (2-6) 
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Pee tecomatac Control 

The automatic control system used to maintain each ship 
in position alongside is assumed capable of deflecting the 
rudder and changing propeller speed in response to measurement 
Peeenceonip s heading angie error from equilibrium, Ay, and 

. 

the separation parameter errors from equilibrium, Aa and AB, 
as well as the rates of change of these variables. The rudder 


control is sensitive to v, wv, B and 8; the propeller control 


moO Cc” and O. 


a) 6 k, Avtk, t+ K3AB8 +k, 6 


kK pAatk a (2-7) 


yeaa 
(where Ky are the control system constants) 
There will be time lags, (At), in both control systems, 


representing the lapse between ordering a control response 


and attainment of the response. Thus rudder deflection at 


(wimtere t = £ +A). 


time t 1S proportional to wy at time t 1 


ite 
Therefore: 


S(t) = kW (t-At) +k (E-At) kB (t-At) +k, B (t-At) ae) 


In the linearized equations of motion w and 8 are 
functions of t. Thus these variables in (2-8) must be 
obtained as functions of t, which is done through another 
PevylersexOansiOn, Carried out in detail in [A~3]. If the 
time lags are assumed small, only linear terms need be kept, 
and the automatic control system equations for 6 and n 


become: 
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oD 
oo [k, +(k.-k, At, )0 “k At, B JAW 


2 
* (oe, k,At,) 0 k At, 0 “JA8 
7 5 E 2 7 
by) m= [ko +(k.-k,At,)0 k At i }Aa ( 2-9) 
(Where Neweme ©peracer o/ot; e.g. rie = 97/dt?) 
Here At, is the time lag associated with rudder response 
and At. that associated with propeller response. 


These expressions for 6 and n are substituted in the 


equations of motion. 


Vee 2oucartions Relating a and 6 to u, v and vp 
Examination of Figure II-2 shows the following 


relationships to hold: 


a) a 


tl 


= I ‘ , — , 
u,cosy, u,cosy,+Vv,siny, VpSiny, 


b) 8 


tl 


u,Siny,,-u,sinp,+v,cosy,-V, Cosy, (2-10) 


Since vy and ve are small, cosy, = 1.0 and cos), oe oe Org 


ince << ee il 
And since Vn SSUQ and VpssuUge and sinp, and sin, are small, 


Ene products V,Sinv, and v sinv, are second order. Thus 


B 


(2~10), a), becomes: 


as UU, “Uy = ((a,) pt4u,)- (Cay) ptAu,) 


where (un) 9 a5 (Un) gs and: 


a = Au,~ Au, 





A Me AS 


Also, sinc€ Yn and We are small, Sin), = Wo and 


sin), “Wp: Thws, (2-10), b), becomes: 


ux. 


— 


d2 


8B tAv,.-Av 


Up_Y RUA va A 


s Vp (Ug) 97a (Ug) (FY pAup-Y, Au, +Av,-Av, 


Zana OGde»r 


Hinererore, Ene equations, (2-10), become: 


a) a ~Au,,-Auy 


b) B zu Ap,-u Ap, tAv.-Av, (Zoe) 
\onrny «/ 
may} 


RE ER I ER ST RE RT EN 


iwo Ship System 


Incorporating all the foregoing discussion, the linear 


eguations of motion for the two-ship system emerge: 


Ship A: | (2-12) 


+ : 2 “L A e 
a) XO AatxX. AB x Au, +X, U, Xx Av, +X. V 


+X* 
A Ba A A A re 


.; 2 
wemmeee ee (ke +(k, -k. At. jf) -k. At. JJ“1]AW 
oeee ret, 2, 2, Lf ee, A 


6 a 3 


+X, [k, +(k, “k, At, )O-k, at, J 21 AB 
a rn ‘nxn 3a ta A ta 


6 D 2 Ee 


HX [ke H(ke kg At, )O-k, at, O 7Jao-ma, = 0 
A A A A A A A 


Se 





MN 
[ 


tL 


> a 
| : 
b) Y Ao-rY, ABt+Y Au .t+¥. U.+¥ Av t+¥. vty: 
On Uy Uy A Uy A Va A Vn A Va A 
2 
+¥. ,+¥. {k. +(k, -k, At. Jf -k. At, fi “JAy 
A 
7 on Ty 2g dy dy ao = 
2 
+¥. [k, +(k, -k. At. )0 -k, At, f) “JAB 
3 
on an Ag 3Q dy on 


ae [ky HK —k 


= At Reece it: laa 
A A A A = A A 


-m, {V,+u_V,4+x_ §.]=0 
peek 
OTA Tn A 


—_ 
c) N. AotN, ABAN Au,+N. OG |+N_ Av_4+N. VIN: qf 
On (Pa Uy A Uy A Va A Va A ba A 


—— 


u 1 


Blais ky #(k5 ok, At )D Bs, ee 5) “100, 
A ei Poa A -A 


Wy [ky +(k, -k, At, )O-k, at, & “lag 
A 


On a A 3n ds A 1, 
2 

ote At) ek. At fF “|Aq 

eee eR eR OR, 
oe pm AG (v, tu) =0 

7 

a) 
2), Precisely Similar Equations for Ship B. 
£) 
Gera. = Au, \u, = Ba 
h ; = WUD ~ -- _— = 

ee ae oY Av. Av, Se 


38 





Meideneitying Zero-Value Derivatives 

Pawan eacomecyumetty Of the Ship about its centerline 
plane, Abkowitz demonstrates that Xs the change in 
x-direction forces of the shiv caused by a change in v, 


with all other variables fixed at their equilibrium values, 


must have a value of zero. Similarly, it is demonstrated 
Eiat Xx = i = ao reese | eos t- 28 ang I-30)... Due to 
the port and starboard symmetry of the rudder, X 1s also 


Zero. Geode og. 1-67). 

Further arguments from symmetry demonstrate that in the 
Seen water Case, Ye Yas N, and Ne aie walieZerO, - oul Ot 
the two-ship system of interest, the interaction side force 
and yaw moment are affected by changes in speed. Thus a and 
N, cannot be assumed to be zero. The effect on this force 
Sao moment, Of a Small acceleration, is, however, still assumed 
negligible, and Y. and Ne are assumed £oO be zero. 

On examining the hydrodynamic derivatives due to the 
orientation of the ships with respect to each other, further 
identification of negligible derivatives can be made. The 


derivatives xy and X representing, respectively, the 


g! 
change in x-direction force caused by small changes in 
longitudinal and lateral separations of the ships, will, at 
most, be small. Experiments have shown that interaction 
phenomena, while affecting Y forces and N moments, have no 
measurable effect on the forward motion of the ships, or on 
ieee Gecistance forces. Call iA-Z|) Thus, XO and X, are 


B 


pemmegual tO Zero. 


Bo) 





te) 


Summarizing: 


V W 6 a 6 u u 


With this information, the equations of motion are simplified 


X =X. = Xt = X% = Xe = X =| KX, = Y. HF Nw = YO (Za? 
oo “Y 


Ship A: (2=14) 
a) X Au,+X. &.+X [k. +(k. -k. At. )§ 
un A Uy A ny Dy On a a 


2 ; 
See Nite DH JAc=m_ au. = 0 
6n 2n RA 


soe ' - i 
b){Y AgtY, AB+Y Au ,+¥_ Avi+¥. V4Ye vo 
On Ba u, A Vx AV) A va A 


me 


> 
+¥. ity, {k. +(k. -k. At. )ff -k, At. & “Jay 
Vn A on 1, 2n 1, 1, aN Ly A 


WV, (k, +(k, “ky Mt, WB -k, At, O 1 A8 
A ol 


+Y, Uk, tke “kp dt, )M-k, At, fet 
A A TON A 
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c) Ny ee Au, tN, Av tN. VatNy Vy 


A A A A A 


| 2 
+N.. +N, [k. +(k, -k, At, )ff -k. At, I) “]Ayp 
YA 6, 1 oreelnts <1 Cae he A 


A 


2 
+N, [k, +(k, ~-k. At, )§f -k, At. £) “]AB 
on 3a An on In An tA 
WN [ke +(k, -k, At, JL -k, At, B yA 
eA he aD, 
=i Ya7m,X (7 tu vy) = 
A JA 
qd) 
e) Peceuecly otintlar Equations for Ship B. 
f) 
g) Ja = Au, ~ Au, 
h) J eg = uoAv, = Ui Avy i Av, = Avy 
ron hea ranging, Equations, of Motion 


At this point, terms are gathered so as to have each 
variable multiplied by the sum of all its coeffieients in 
each equation. 


Remembering that (variable) = @(variable)/dt, then: 





Buu) = d(Au) _ d(u-u) _ ou O. (2-15) 


one 
But Uy PoMeOms ean... SO Hau) = OGL 


Al 





Lives: 
Sigel” 2. (2— 16) 


2 
Pipe eet kee ek. Ac. Nf) -x «6k. At. U “jAq 
eon TN One HK Ny, 6, 4p 


+ [X. +(x. -m,)fJ JAu. = 0 
Uy Un A A 


2 
Pelee eke yey (Kk. kK. Awe YS -y k. At. B “JA 
A ™ -~ MM Sn A 2a Ny Oy 2y 


2 
Coe oe GO BY k At AB 
ar ne a Al A ie Ch) ear 


V 


TR 9," yethtgue™, 2 Av, 


+[Y. k, +( (Ys -m,u_)+¥, (kj, -k, At, ))f 
On 15, Va AO On 2n ly 1 


2 = 
Pee ar St ar Ms JAp, = 0 


2 
Smee tN Kk. )4+N. (kegeloedt, )BieN bk. At, J ~ ldo 
Oy Ny Sp Ng On Mp 2Q Ny Sn 4n 


4 3 il 


+[(Ng +Ng kz )4Ng (ky -kz At, IO -Ny k, at, 6 21 AB 
A A A A A A A A A A 


+N Au_+[N_ +(N. -m_x. )0 JAv_+[(N. k. +((N* -m_x_ u_) 
Uy A Vn Vn A Tn A On ly Wn A 5, © 


See 2 1 


eek. -k. At. ))D +((N. 
on 2a ae Vn 2y A “A 


, See ee ae “Jay, = 0 
A A i A 
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In writing the equations for Ship B use is made of the 
relationships provided by equation Geena (2-14), to 


substitute for Au, G@eevemin the whip B equations. From 


B 


these equations, Au, = Au, +fJo and 


Av, — Le +f Av, “ UoAbn = UAE: dhe bise 


Simson: (216) 


aja. ee (xe, =I At. Jax. yO 
Mp yp Np SR Mp 2—R Up 


+ ((X, “M,) x) k 


Z 
¢ At, en) ae 
B B £6 B 


WRX, my JAu, a0) 


rime mw wy. kK. )+(Y,. (k, -k, At. )+Y._ )& 
BR 8p 3p Sn 4—3 33 lp Vz 


2 : 
Ges ome) -Y. k, Me _ Peieraie | Au, ay -m,,) 0 JAv, 


B on 4_R ot B B B 


si lugty to (Ye mye JAW, 


2 
—u (Y. - (nee _ k t = 0 
UG v mB ( ( bp MEX gy) ae a JAY, 
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if ne zs 
Per . rel Ke )+(N (kK 6 ky 


B pe b B B B “B 


tN kK, ) + (N, (ky “k, 


+LIN, 
B B & B B B B 


2 
-N. k, At. )){“]A6+IN 
Sp 4p 1, aS B 


Nemec (N. 
ony @ a 


-m_x )f JAW _+[[N, k 
B Be eB A" 85 


i 


creciewaitex WU )tN. (ke —-k_ At...)-u 
Vp B Jp oO On 25 1, 1, O 
Z 
Cpe )-Neek. At. Oi -lAta= 0 
vp Zs 6. 25 1, B 


I —N 

At. )+n )O . k 
+-N 

At, ) +N. 


JAu,+ (IN, +(N. 


gies 2) Ao 
B ee 


9 +((N. -mx_ ) 
B 3 


“UN, ) 
B 


These, then, are the six equations of motion in six 


vwerriables, 


are of the form: 


(Axl) Aa + 0 =: (AX3) Au, AP 0 

(eri )yAa + {AY2)AB + (AY3) Au, + (AY4) Av, 

(AN1) Aa - (AN2)AB + (AN3) Au, 7 (AN4) Av, 

(BX1) Aa + 0 + (BX3) Au, F 0 

(ey )Ao + (BY2)AB + (BY3) Au, a5 (BY4) Av, 

(BN1) Aa + (BN2)AB + (BN3) Au, + (BN4) Av, 
Mio Sea Set Of Six Simultaneous, 


Dimes? <n nowns. 


equations 


Weea tO denote: 


iMeMeott taney storm fOr SOLutilvon . 


ne term BxX3 , 


The equations 


(2-17) 
+ 0 + 0 =0 
+ (AYS) AY, + 0 =() 
+ (AN5)A), + 0 =0 
+ 0 + @) =0 
+ (BYS) Ap, + a aaa 
+ (BN5)Ap, + (BN6) Ap, =0 


linear differential 


for instance, 1s 


Paros a Kw emuataon, 3rd term, (or Aun term) . 


Meer the coefficients is, itself, a polynomial in DD aes 


44 





Prampemceen trom equations (2-16), for which (2-17) is a 


Shorenand notation. 


Pte 2Olving the Equations of Motion 

Ef the terms AX1, BX1, etc., were just numerical 
coefficients, straight-forward algebraic technigues could 
be used to solve this system of equations. Yet, while these 
coefficients are, themselves, polynominals in the operator J ; 
PMieovreZ (Shows that J) can be treated as an algebraic 
Calantity . Gites, pew L-26,) 2/)> Thais 18 possible because 
the hydrodynamic derivatives are defined as the slope of a 
Beees Ot MoMent versus a dynamic variable taken at the 
Selenium Condition. Thus the terms in the coefficients 
Peter ecrte,, Other than RG, , are constants in time. 

Wiewsotucion for any of the variables is, then, the 
determinant of the coefficients, with the column corresponding 
to the variable being solved for replaced by the right-hand 
side column vector, divided by the determinant of all the 


Goeetricients. 


fOmee x<anp Le: 


(See next page) 
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OO -AxS~ 0 0 0 
OPA gery Ss AY4 AYS 0 
0 AN2 AN3 AN4 AN5 0 
OOD eas —0 0 0 
O BY sb ioe Pye by > BY6 
O BN2 BN3 BN4 BNS BNo 


Aa = | ~ TDET | 


AX 0 AX3 _.0 0 0 
Peeiv 2 AY S AY4 AY5 0 
ANLAN2 AN3 AN4 ANS O 


Ext 0 BxXs? 0 0 0 


(2-135 
Ex BY Bye pr4s BY5S BY6 
BN1 BN2 BN3 BN4 BNS5 BN6 
Bvaemation Of DET leads to a polynominal in ead 
me Oo, Sums Of products of polynominals in i , SER nase ro; 
etc.), whose coeffieients are various combination of the 
hydrodynamic derivatives and control system constants. Thus 
the solution for each variable would be: 
: - 0 
variable = pe ee eee 
a, ats a. D Oe cit (Oa195 


[ie "beGero Other than zero, the solution for all variables 
would be identically zero for all time after a small disturbance 
Paegeciulitiocium—---a physical impossibility. Thus, the 
denominator, the determinant DET, must be set equal to zero. 


Setting this polynominal in 0 , resuiting from the evaluation 
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Poaioeoowel torzero, the roots Of the polynomial can be 


obtained, where Oy1 OO Bia C are these roots. Then: 
Wwatieanale = | | SO 
(DB ae) is’ ieee) Beene (# Ga) (2-20) 


Abkowitz demonstrates how the operation 7525) [21 
is equivalent to: 


THz) (2! = eo fe “fay at (cao 


0 Ons =O) te = 
Thus +75 =) =e fe (O)dt = c,e 


OE 


Carrying out this operation successively to evaluate 


equation (2-20) yields: 


Vaenaole = ¢c eae + C Bo Sea ieee ert cen 
il zy n 


~_— ow = 


where the c, are arbitrary constants of integration. 


Pee seco tiching Che Stability Criterion for the System 
Since each of the variables is a sum of exponential 
fermms an tame, it es clear that the only way the variables 
can all go to zero as time increases, is for all the O; (fons 
all real parts where the O; are complex) to be negative. The 
criterion of stability of the two ship system is for all the 
Variables, representing, as they do, deviations from the 
equilibrium condition, to go to zero as time increases. For 
when all the variables have gone to zero, the system has 
regained the equilibrium condition, meeting the definition 


of a stable system. 


47 





meee stability Gf the system requires that the real 
parts of all vg; be negative, As the ¢; are the roots of the 
Sieeacteristic equation of stability resulting from setting 
the polynomial in Tf , (which in turn results from evaluating 
DET), equal to zero, the stability requirement for the two-ship 
system is that the real parts of all the roots of DET be less 


Eilean Zero. 





Caapcer FTI. 


Peececaure 


ie Background o 


Pemetei-ea at tae conelusion of Chapter I, a modified 
OLNA/MARINER class merchant hull is the ship used in the 
Mathematical model. As indicated there, this was done because 
open water hydrodynamic derivatives, in non-dimensional form, 
were available for a MARINER, ([A-4], [A-5]), and because the 
only good, quantitative interaction phenomena data available 


Peaentor a Ship Quite similar to a MARINER---the R.F.A. OLWNA. 


Mee oben Water Hydrodynamic Derivatives 
Reference A-5 presents experimentally determined values 
for most open water derivatives, measured with a planar motion 
mechanism, for a MARINER model. Reference A-4 utilized these 
Values, with the exception of Xo for Which a ditierene value 
was determined. the accuracy of the values used in [A-4] was 
borne out with full scale sea trials and these values are used 


in this thesis for the open water derivatives, for ship A, and 


wee lacreea in Table ITTi-l. 


Hiz. Obtaining the Interaction Hydrodynamic Derivatives 


N,N 


ae Yu! Yer zs g° 





Values for hydrodynamic derivatives, as such, arising 
from interaction effects, could not be obtained. The inter- 
eemtonmedata for OLNA consisted of dimensionless forces and 


Mem@emes at Side-to-side ship separations of 50 and 100 feet, 
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for all longitudinal positions of the midships sections of 
the ships between -600 and +600 feet. Also available were 
force and moment data for OLNA at side-to-side ship separations 
Of from 30 to 150 feet, with the ships abreast. [A-1] (See 
heegqurces —I-3, 1-4 and I-5) 

Using the data illustrated by Figures I-3, 4 and 5, the 
g! Ny and Ne could all be 


Getermined were at side-to-side distances of 50 and 100 feet, 


only positions at which Yue x 


with the ships abreast. (Feommis 1 and 6, Figure ILE-L). Gage 
limitation of the examination Of System stability to -enly two 
ship positions was felt to be too restrictive. As the only 
means available for obtaining the desired derivatives at other 
positions of interest was extrapolation and interpolation of 
the OLNA data, this course was followed for both Y force and 


-—- 


N moment. As data was available at side-to-side distances 
food, tor Short), of 50 and 100 feet, these values were chosen 
as bounds to the range of interest, just as values of ~-600 and 
+600 feet were chosen as bounds in the longitudinal direction. 
The way in which Y force was interpolated in the range of 
interest will be illustrated. 

The Y force was known along three lines of position, as 
Sem@ewnein Pigures 1-3, and I-5, and as in Figure III~-l. To 
determine Y at positions other than those known in Figure III-l1, 
it was assumed that the ratio of change in Y between 1' and 2' 


to the change between 1' and 6' was the same as the ratio of 


the change in Y between 1 and 2 to the change between 1 and 6. 
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Pee ee a re) (YY) 

eee oe Wee eT at 

imoommeo i etnzougin 6 im Figure Til-1 were chosen at 
[MomroOot ameervals, and the interpolation for these values of 
STSD=60, 70, 80 and 90 feet was performed at 50 foot intervals 
Mime tame, vyalues of Y were obtained at all intersections of 


a grid with lateral spacing of 10 feet and longitudinal 


spacing of 50 feet, as shown below: 


eal eo SS has 50, 
| 
ee | = =.$7S pe ee 
{ | 
| § 
—-4dt-— pe ae 
' 


Jal copiers) BAL Aba 


teemeesills Of Such interpolations fOr both Y nd N are 
shown as figures III-3 and JII-4, where the curves for STSD = 
50 and 100 feet are from [A-1] and the others are interpolated. 
The interpolated curves are felt to give a reasonable picture 
Of the way the interaction force and moment would be expected 
to vary as STSD is changed. And, of course, the points on the 


line a=0 are the same points as those plotted in Figure I-5. 
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Obtaining the interaction forces and moments.over this 
meme Of 50 oeStSD < 100" and -G00° < a < +600° is, though, 
Beecmarstep alongutne way tO Obtaining the derivatives Yue Yer 


Ny and N which are what is really desired. 


B! 

The longitudinal separation parameter, a, is represented 
(ireeGgmely On the curves @©@f Y and N. The lateral separation 
parameter, §, is related to the side-to-side distances shown 
on the figures by a constant, for a given pair of ships. 


Since $6 1s the lateral Separation of the centerlines of the 


Snips, it 1S apparent that: 


8 = STSD + (BEAM + BEAM 


Ship A Shi ee 


i 


erss + Constant (32) 


HiMcpeira: sewO Sihies Of 50" beam, an STSD of 50° is 
se valeme to B=100°. Therefore, Figure III-3, for example, 
represents a surface in space, (as sketched in Figure III-1), 
whose horizontal coordinates are a and 8 and whose vertical 
coordinate is the value of Y everywhere on the surface. Then, 


E@mget, oY faa, Or Yue at any point, the slope of the Y surface 


Meret omardmrectilon, at the point, must be obtained. Similarly, 
Ye is the slope in the 8 direction. For example: 
ay tome e— Y CL00,60 ) 


((¥ (150,60)-¥ (50,60) ) / (150-50) 


co) ae NC, 2) NGO 0760)) 


ee eeeeer 0) —-Y (L00mS0)) (70-50) (ee) 
im this way, Y, everywhere in the range (50'<STSD<100'; 


55 





Pa wom 50 ), Was Obtained, as was Y, everywhere in the 


6 


memoem cod ~olob<90°; —G600"= a < 600°). (The values of yy and 


Yp at the extremities of the range of known Y values, in the 
Peirce Ten OL Citrerentiatcion, Could not, of course, be obtained 
by this simple procedure). 
| The values of Ny and Ne over the same ranges were obtained 

in the same way. 

A computer program, program INTCOF, was written to take 
the values of Y and N at the locations known from Figures I-3, 
#4 and 5, to interpolate to get Y and N throughout the range, 
ememtchen determine a Yas Ny and Ne in this range. The curves 
Pm Ljunecmitm-3 and Lii-4 for side-to-side distances of 66, 
(ye oO ane 90 feet are plots of the Y and N outputs of the 
INTCOF interpolations. 

Stee chee cURVes in Figures I1I-5 ana ILIi-6 for 
abe o0, 07a) and Umfect are INTCOF outputs. And in 
Figures III-7 and III-8, the curves are all INTCOF outputs. 

Program INTCOF is described in detail in Appendix II, 


and the results of the program for the modified OLNA/MARINER 


ship are included therein. 


Bae <7 Nae 
u u 





AS newed in Chapter 1, e@sperimental investigators in the 
field of ship interactions have noted that interaction 
phenomena vary with speed in the same way as the resistance of 
the ship---or that interaction effects vary as the square of 


the speed. The OLNA data used for interaction effects was 
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poLomwmea at a Speea ot PS 9knots; which was also used as the 
equilibrium speed for the results presented in this thesis. 
(However, all force, moment and derivative valucs are 
presented in non-dimensional form and are Toe at other 
Speeds.) 

Remembering that Ye for example, is dY/du with all other 
variables at equilibrium values, and knowing Y for a given 
Position at 15 kn@ts, it is a Simple matter to get Y at 14 and 


h6 knets , and then = at. Ism@iene@ts:. 


meets, = v(ts) & (16/15) - 


i 


mye = 15) & 14/15)" 


e) ¥,, (45) =v ro) ara) ) 7 CG -14) x. 6378) ) (oy 
Ny is obtained by a similar process. 


GF Xx 7 Y I N e 
n n Ty 


mee ee 








The effects on forces and moments of a change in propeller 
Speed are also required, and are represented by xl Yn! and Nn 
The derivative x is determined by using the known values 


of Effective Horsepower, (EHP), and RPM for a MARINER at 


speeds around 15 knots. 


a Wes) = (BHP tr) x550) /(15x1-6878) 
i) Saenoy = (EHP (16)55500/ (16x1<6878) 


c) x (15) = (X(16)-X(14)) KA( (RPM (16)-RPM(14) ) /60) (3-5) 


ie eavelie CL x So obtained is in dimensional form, and 


1s then non-dimensionalized to permit broader applicability. 
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Piemomnele serev lABINER requires a 1.2° port rudder 
angle to maintain a steady course in open water. {[A-4] Since 


mieevalue OL Y. Lor the MARINER is known from reference A-4, 


6 
it is possible to determine what side force is being generated 
by this 1.2° rudder angle, which is denoted by Uo eee When 
Maintaining a straight course in open water, this rudder side 
force is offsetting the side force caused by the asymmetry of 
a Single screw. Thus, on a steady course, Geen 

The dimensionless Y Given in [A-4] can be converted to 
Peeeemsional form, Which iS in units of lbs/radian. The 
actual Y force exerted can then be calculated by multiplying 
the dimensional Y by oP ol Sivindsasside Lorcemin 
pounds which is the negative of the force being caused by the 
propeller. The force exerted by the propeller can be expected 
to vary as the propeller Speed squared, and, thus, ions etane 


propeller side force at a given speed, it is possible to 


determine how this force varies as propeller speed, n, varies. 


a) Sara = Y ¢XNDF (@15 kts) x (op eed radian) 


Ib) ¥ (hSekts) = ~Y 


PLEOD rudder ‘+7? 


uy (ole) (ies) x ALA) FAs) ec 


oe 
Een DEO 


Ga (16) Y 
prop PEOD 


(15) x (n(16)/n(15))* 


e) vy, (45) = ron ae ee ee 


te 
(Where NDF is the non-dimensionalizing factor for Yo) 
The value of es is similarly obtained, and the YO and 


N, So gotten are suitably non-dimensionalized. 
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A computer program, program DETERM, was Wiese tenet 
Getermine the values of vy! Noe Xe aa and No: The inputs 
for the program are the values of Y and N at the positions of 
interest, gotten from program INTCOF, and the EHP, RPM, speed 
and BO eiecs values for the ship being examined. The values of 


Xs and Ns in open water are also inputs. Program DETERM is 


described in detail in Appendix III. 


Reemeecalciudmation Of Interaction Derivatives for Ship 8B. 

All discussion concerning interaction hydrodynamic 
derivatives, so far, has been for the OLNA/MARINER ship which 
is assumed to be Ship A of the two ship system. Ship A serves 
as the reference for measuring both longitudinal and lateral 
separations, o and §, and is assumed to be the overtaken ship, 


MOemred LO the jeft of Ship B, the overtaking ship. 


Fae a era Sa ee hr 
/ ~S B, 
\ B N—--- - > +> 
N ae 
ee 


“ee ee 


This coincides with the case for the OLNA experiments and all 
Blots of the interaction variables vs. o are for the overtaken 
ship on the left. 

But hydrodynamic derivatives for both ships, at any 
relative position of interest, are needed to write the equa- 


tions of motion for the system. Since the modified OLNA/ 


Ge 





MARINER used as Ship A is the only one for which it was possible 
to obtain all required derivatives, Ship B must also be 
Ponsidered to be this type of ship, in order to be able to 
Obtain derivatives for it. If the two ships are identical, 

Mae £Orce Or moment felt by Ship A at a given a and 8 are the 
Megataves Of the force or moment felt by Ship B® at the same 
peeaceee-c%, TO illustrate: 


ie 2. 





—— o— em eee oem oem oe eee — 


PeOmmonmewseac 2 comalis -Y on Shido A at 1. If Ship A is 
Betractea to Bin position 1, then B will be attracted to A 
in the same way at 2. But because of the reference system, 
Tieewetosecatea i ragure Il=2), an attraction for Ship A’ is 
Meeesltive Y Lonee, while for Ship B an attraction is a 
negative Y force. 

Thus, in program DETERM, the Y and N values entered for 
Ship B are the negatives of those for Ship A when a=-a. As 
an Gxample, with STSD = 70' and a = 300', the non-dimensional 
Moevotrcump A ws 1652 x ioe The Y value entered for Ship B 
iiwindas case as -.1689 x TWO which would be the value for 


Ship A when a = -300'. (These values are from the INTCOF 
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Meepucwodui or Appendix If.) Wath the inputs SO adjusted, 
program DETERM will calculate correct values of yy and ae £or 
moth ships. (ier calculearrons for Xe ve and No are indepen- 
fene Of thse considerations). The output of program DETERM 
me, cthererore, AW, ¥, and No for Ship A, (which are assumed 
Eo Dewthe same for Ship B), Bs and N., for Ship A and yy and N. 
Poe, 1G it were in Ship B's position. 

To bypass the necessity for considering both ships 
identical, provision was made to let Ship B be of a length 
different from Ship A, but of the same form. All form 
coefficients, plus the length/beam and beam/draft ratios, are 
assumed the same for both ships. This assumption implies 
that the non-dimensional values of the hydrodynamic derivatives 
Bom yonip Gb aveuthe same as those for Ship A, (with some sign 
pee ence... as in Ye and Ne), when all derivatives are 2 
dimensionalized with respect to the length of the ship to 
which they belong. But for the two ship system, the charact- 
eristic length used in the non-dimensionalizing process must 
be the same for all derivatives. It was decided to use the 
length betwoen perpendiculars, LBP, of Ship Ps Rhus jo. 2 7 
Bens Or a Tene ch G@poeeomeme chan that of Ship A, the dimension- 
HesseGeriVilives of Ship A must be multiplied by (LBP),/LBP ,) 
to get the viimensionless derivative for B with LBP, replacing 
LBP, in the non-dimensionalizing factor. 

Tiemiieditieation Of the Ship A derivatives to fit Ship B 
1s done in the early part of a third computer program, SOLVE. 


(This progiaum is used to actually solve the equations of 








BeErem for etn Gtability FOoOts, as outlined below, but modifies 
Bile"CeriVatives, ac Stated, before solving the equations). 


SOLVE is described in Appendix Ivy. 


V. The Complete List of the Nydrodynamic Derivatives 

coe til lists the de@ivatives from reference A-4 for 
SteeereiNER and these are Che™values used for the open water 
pelivaeaves for Ship A, and whitch are modified to apply to 
phaip B in program SOLVE. 

ijl eii-2 lists the demivatives associated with propeller 
speed which were calculated by program DETERM for Ship A. In 
program SOLVE, they are used for Ship A and modified for Ship B. 

The values of Yue Yee Nv Nes vy and N, are position- 
a Popeeeewecence the values Giy and N acting asa result of 
interaction effects vary as thr ships change positions relative 
memeecmocher. in addition, @Wltese derivatives, “even™it ships 
Pema sare identical, will nol be the same, for the two ships. 
Seolveoertl-3 gives Values of tlisse derivatives for the two 
Phips, assuming them identical, for two typical positions of 


the ships relative to each other, 
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Sie LLi- 1 


Derivative Non-Dimension- Non-Dimensional 
(Selizang Factor ...... . Value 
(X,~m) W2eL* ~840.x10 > 
2 ue =—2 
x WOE, 120.x10 
y 1/2pL2u -1160.4x10°> 
V O 
(Y-m) 1/2pL3 -1546.x10"> 
ae 3 _ => 
(Y¥emu_) VW 20k Uo 499.x10 
a 20, 8. 6xlome 
Y, 1/2pL?u2 277.9x10 > 
N Oot a Cees 
V O 
ce Wa ole De Tei 
— 4 =— e 
Nyx ou.) 1/2oeL ve VeGex LO) 
(Ne I.) 1/2pL° penis 
N, iy 2a -138. 82108" 


mass density of sea water 1.9905 lb-sec’/ft’' 


ie 
i 


Pigeons A, Li ft 


u = equilibrium speed in ft/sec 


TABLE III-2_ 





Derivative Non-Dim. Factor Non-Dim. Value 


1/2pL?u, MG on 
1/2pL7u, = 52alone 
i =o 
1/2pL Vo a2 60 
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ee a ce 


ie Chapter Ii, 


lie Wee 


Taba LLi-3 





o=0' i STSD=60' 
Deriv. Nien am. Fackor Ship A 
se Non-Pam. Value 
Y 1/2pL2u Nepe Selo) > 
u O 
~~ 3 _¢- aS 
ae L/D Uo yell 
: 2 =o 
XO Ive ols ue - 16x10 
2 a aS 
i 7 Zou us .87x10 
pieho =. 
v 172i Us ein 0) 
ee = 
ve I Foe ue ~42x10 
Paw 400'; STSD=80' 
e — = Gone 
N,, Miesxi0 = 
= -_ 
is, ) as .044x107> 
oN“ Ps a -5 
Y : 
g A oe 447x10 
N ~.068x10> 
Ne Een? 
Pee lerned of Solving for Stability Roots 





Sidabjon de 
Non-Dim. Value 


ee 
52550000 


G7 eon 


igen: 


.87x10 > 


5 


.126x10 





-16.7x10~ 
~.093x10> 
-.111x10— 
-~.035x10- 


1s 
.098x10 ~ 


oe aan: a 


shown that the stability Cid Beem 


teristic equation for the system was obtained by evaluating 


the determinant formed by the coefficients of the variabies. 


These coefficients are polynominals inD, and it was stated 


that the operator could be treated as an algebraic quantity, 
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so that the result of the determinant evaluation was another 
polynominal inJ, whose roots were the stability roots. 

The coefficient determinant to be evoluadeeduis a 6 xX 6 
determinant, given as equation (2-17). The determinant is 


evaluated by pivotal condensation (A-13, pg. 155, problem 


4-21). An illustration of pivotal condensation follows: 
aay aio Ue ae Ain i ay ee in 
a a 
Ao ao 5 Soe toss ao ae flea 
= Malco Cleo. 6s 
ihiewie 21. 23 arn 
ani a5 eye cy! aoe ant Anos a 
iL Oa nae 
! ! } t 
Anz Aggreces aon Angrrre: an, - 
=A. = 2 ! ' 
ae | : a: oe a (3-7) 
a Ei ie Ghee: a 
la az nn 


Where a oa a,.~(a,./a,,) (a,4) 


The procedure would then be repeated with las+...... a 


The matrix properties whereby the interchanging of two 
rows or two columns changes the sign of the determinant are 
also utilized. After three applications of pivotal conden-~ 


Sation and use of these two matrix properties, the determinant 
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in (2-17) becomes: 
(AX1) (BX3') (BN6)*|AY4 AY5 AyY2 


Ba By BY 2" 


AN4 AN5 AN2 (3-8) 
Wirewe +) 6x3. = Bx3 ~ Bx (AX3/AX1) 
BYA = byt — BY6 (BN4/7BNG) 
BYS' = BY5 -— BY6(BN4/BN6) 

BY2' = BY2 - BY6(BN4/BN6) (3-9) 


Sain rngsene Gvalluation of (3-8) to ats conclusion, and 


Metis zing the relationships of (3-9), the determinant becomes: 


Dit = 1 (AX1) (BX3)~-(BX1) (AX3) ] 

x{[(BN6) (BY5)- (BY6) (BN5)]x [ (AY4) (AN2) - (AN4) (AY2) ] 

+ [(BN6) (BY2)-(BY6) (BN2)]x [ (AY5) (AN4)-(AN5) (AY4) ] 

+ [(BN6) (BY4)-(BY6) (BN4) 1x [ (AY2) (AN5) - (AN2) (AY5)] } 

: (e=ng) 

Where the various terms are defined in Table III-4. 
Program SOLVE, after appropriately modifying all hydro- 
dynamic derivatives to apply to Ship B as well as Ship A, then 
falls a Subroutine, DETER, to Garry out the mathematical 
operations to evaluate (3-10) at one or more positions and 
with one or more sets of control system constants for the two 
ships at the various posterons,. Bach Of the terms an (3-10), 
ee, BNG, AY2Z, etc., is a polynomial in ) , and the various 
multiplications, subtractions and additions are carried out in 


seorouei ine DETER by using IBM provided Scientific Subroutine 
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Pacwege (Sor) sibroutanesior arithmetic operations with poly- 
Momials. Program SOLVE, (and all the subroutines it calls), 


is described in detail in Appendix IV. 


ere serene mee ee eS rm rr ee 


Reference A-l1 discusses automatic ship control and 
utilizes some typical values for the control constants in a 
system sensitive to wW and w. Sensltigity LoOvansesror 
Pe oemerer, Such aS YW, 1S proportional control, and to an 
Seer Darameter rate, such as v, PS ratemeoitro]. “lve 
equation for rudder control in such a system, (not 


Somsidering Signs); is; 


6 = k)p tkyy (35 )15) 


ioe jA-l1|, Schiff and Gimprich use values of Ky between 1.0 


and 5.0, as typical values. The only values used for k. are 


Z 
Ormor 1.0. 


With Ky Sled al" neading error will give rise to 1° 


ee correcting rudder, while with ky = 550; a ieee ror yo 


@amse a 5S” rudder correction. If Ky = 1.0, a heading error 


Bate on I-/7sec. gives a rudder angle of 1° for correction. 


Thus, Ky was considered to vary from 1.0 to 5.0, (with smaller 


values investigated to determine some trends), and k. was 


Z 


varied between 0:5 ana’ <4. (Defining 6 and yw to be in radian 


@@eesenot affect Ky and Ko). 


In determining ranges to be used for Ky and Kye 


determine sensitivity of the control system to lateral 


which 


Separation error and its rate of change, 8 and 8, no basis 


12 





mee found 1m) che literature. 


Thus, these constantS were 


estimated from what the author expected to be reasonable 


corrections of the rudder angle for lateral separation 


Sierors. 
Signs are not considered): 
—— kB “+ kB 
Based som the oallhnor s 
it was assumed 


operations, 


Paemoe a rudder correction 


Rudder angle is sensitive to 8 as below, 


(again 


(3-12) 
experience in underway replenishment 
Eidtadn Orcor 4 WO eet an 6 would 


Srenavoue i . Ane a ehange. in B of 


Mae SOc. was felt to call for 3° of rudder correction. 
Therefore: 
Dei 1 radian~i + 
eee sso) sys OO! 
eS i waerran ~ 3 ~ = 
eee i 3° ee feet) 
As a result, k, was investigated between .0005 and .01 


and k, between .0l1 and .l. 


A 


A propeller speed control, 


sensitive to longitudinal 


Separation error and its rate, a and a, was assumed in the 


theory. For this sensitivity, not considering signs: 
n= k,a + ka ales 
To get an idea of the magnitude of Ke and Kes estimates 


were resorted to. 


MOUSGeetlan @ would prompt an n change of 1 RPM, 


ie wasoedssumed that a typical error of 


(about 1/4 


Peeemror a tTIARINER) and a |] £t/sec. value for a would prompt 


a change of 4 RPM. 


Therefore: 


7 4 





~ 1 RPM i aeisbe cs: 


5 10 £t ~ 6ORPM/RPS 600 ee 


~ 4 RPM ip eee = 
mt Gk «Cts CC OO yoke 


So Ke wasmascsumea tO bemoeeveen .0005 and .01 in value, 


and Ke Ipecreveen wOL and .2. 


wages Control Constant Signs 
When either ship develovs a positive error in J or v, a 
positive rudder angle is required to correct for the error. 


Thus Ky and Ko are positive for both ships. 


«> 
Wie 


As 8 or § increase in a POSTELVEe daicecetOn, It 4S decareo 
feomeurieottp & LO the left and Ship A to the right, to 
compensate. Thus as 8 increases, a negative 6 is required 
On Ship A, but a positive 6 is wanted for Ship B. Therefore 


K and Ky are negative for Ship A and positive for Ship B. 


V> 





o> @f 





AS a increases (or ~a decreases in magnitude), it is 
desired to increase Ship A's speed and decrease Ship B's. 
Therefore Ke and Ke are positive for Ship A and negative for 
pitip BB. 

Both rudder and propeller response will have time lags 
associated with them. The rudder, however, can be expected 
Po eueapenc more quickly than preselier speedy Thus rudder 
time lags from 0 to 4 seconds, and propeller time lags from 
He-OeG Seconds are examined. (The desire is to determine the 
Way in which time lags affect stability, and errors in the 
actual values used will not affect trends developed.) 

io summarize, Table IZII-5 lists all the control system 
parameter ranges examined. 

TABLE III-5 r 


Ship A Ship B 


5 ia Ky < ae, 
_UCC> < K, S oe 
non iS Ky 2 ali 
eal < Ke < eS 
-.2 < Ke < =O 1 
0 = At, “ ay 
0 S At. < one 
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Chavter LV 


RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


ee MeL Oouction 

The intent of this thesis is to develop and solve the 
equations of motion for the two ship system, and to obtain 
Pieestability roots for it. The ne noel with the computer 
programs used to perform the calculations, (which would be 
prohibitive by hand), permits determination of whether or 
not the system is stable for a given set of control constants 
and time lags at sane positions. In order to examine the 
accuracy of the solutions and correctness of the method, 
Various sets of control system constants, for various ship 
MOstti1ons, were used as inputs to the computer programs and 
the effect on the stability of the system observed. 7 

It is beyond the scope of this thesis to attempt an 

Optimization of, or even complete examination of, the effects 
of variations in the automatic control parameters upon 
System stability. The purpose is to present the method, with 
only enough calculation to assess the value of the theoretical 
formulation and computer solution. As a result, the various 
combinations of controls and ship positions examined are by 


no means exhaustive---but are representative. 


el. Results 
Dh Glew seca of Control constants were evalwated at each of 


four positions, for ships of different size, with and without 


ie 





time lag, and are presented in Appendix IV as typical output 
from program SOLVE. 

To perform a somewhat more complete examination of system 
Stability, a single position, (0, Boe > FObe= - OO) yin Iwas 
egosen, and a total of 91 combinations of control constants 
and time lags was run. The results of these runs are 
Summarized in Tables IV-1, IV-2, and IV-3. In these tables, 
the number of zero roots, number of negative roots and number 
ea peOcmtive roots obtained for the input control constants 
Peominaicated. In addition, the largest root for each condition 
1s listed, unless the largest root is zero. In this case the 
largest negative root is listed. The effects of changes in 
some of the control constants on the value of the largest 


Moot are shown in Figures IV-1 through IV-5. 


mie. Discussion of Results 
The first part of Table IV-1 summarizes the effect of 


varying k with all other control parameters zero and no 


alee 
time lags. These results are illustrated in Figure IV-1l. 


. k : 
Rarying 1 alone, from a very low value to its maximum, 


n! 
decreases instability for a time, (the value of the largest 
moet 1S positive, but decreasing), but then instability 
increases as iL. increases. When ois and ae are both 
increased mong this same range, the system instability is 
lessened, with stability actually acheived for part of the 
range, but the trend for instability to increase for the 


higher k, values is still present. 


1 
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ii wwe iivestigatea to get the results presented in 
this chapter consisted of two identical ships. This 
permitted a check of the method to see that the solution is 


symmetrical for this symmetrical system. As a check, a run 


was made with a, = 3 and pa = 0, which yielded results 
identical to those when on = 3 and i = U0. “huss he Svoeen 


of equations, and their solution, appear to be correct in at 
least this respect. 

The second and third parts of Table IV-1 summarize the 
results of varying kK, and Ke tor the twowships, witcha iil 
other constants zero. As can be seen from the table, and 
mrom Figure IV-2, neither Ky ve@nic Kes 


Bimeoth Ships, is sufficient to bring about system stability. 


alone, even when present 


As kK is increased, for one or both ships, with all other k's 


zero, system instability is actually increased. And as kK. 
1s increased, the value of the largest root is unaffected, 
though the values of the other roots do change slightly. (The 


symmetry of the solution was again demonstrated for these 


constants.) 


k k k k , : 
Next, 3! 3p) 21 and 2p were all given values in 
their expected ranges, while ae and Se were varied. The 


results are summarized in the fourth part of Table IV-1 and 


Mtciiatea ain Figure IV-3. In these circumstances, 


increasing ie and re reduces system instability until a 


Siroaole Condition is attained. 
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TABLE IV-L 





itary s@L Stabt tity ROOtS Obtained for Various 


Combinations of Control System Constants With No Time Lags 





¢@ = 90'; StSsp = 0'- a a. 
eee ee = 526.5 °; Shim B LBP = 528.58 
Pea sted Ke = 0 





Ky # 0 # <0 i: >0 Largest Stable? 
A Res. Rts. RES. Re a 
0 | 0 B 6 2 122710"? no 
.001] 0 1 Zi 2 1201x1077 no 
aoa! |) a0 1 7 2 9343x1077 no 
a 0 1 7 2 3293x107 no 
1. 0 1 7 2 6980x107 no 
oe 0 1 7 2 .7186x107> no 
ay 0 1 7 2 oe 0n no 
A, 0 1 ie 2 .7289x10 > no 
5. 0 1 7 2 TetOxtO = no 
0 ah 1 7 2 .7255x10 °° no 
.001 | .001 | 1 7 2 iy sso) 4 no 
mL) 0d 1 a 2 "G8Gtxton- no 
al a |e a 9 0 eo ie ‘yes 
1 1 1 9 0 teh tO = yes 
2 2 1 9 0 ~.9270x107* yes 
3 S a 9 0 = enna yes 
4 4 1 2 0 ASB aR tO) yes 
5r 5 1 a 0 Se Se elO yes _ 





eee 
nl — SS OS A SS Ee ——— ee 








Oe 


2005 


- 008 


201 


50) 3: 


-05 


Oy 


+ O 
Rise 


1 eee OV =). 


ee et ee ee 

















(Gorm. ) 

i 

one a. eae Stable? 
6 2 .2809x10~ no 
6 2 3720210 © no 
6 2 .5848x10 + no 
6 2 7242x109 > no 
6 2 .8348x10°* no 
6 2 3720x1060 - no 
6 2 3720x102 no 
6 2 .4509x10 7 no 
6 2 .4940x10 + no 
6 2 enon a no 
6 2 .9741x10°? 


6 2 See 
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ieee ay 1 


eee 








(con tv) 

sy eo. | BS] 2] terse | seaptes 
0 2 mip yo - no 
0 2 i no 
0 2 : no 
0 2 " no 
0 Z 7" no 
0 Z ‘ no 
eae ON. 2 : no 
BOO. Z : no 
.008 2 " no 
ee 2 - no 
03 2 y no 
205 2 : no 
Ol 7 2 i no 

3, = -.005 = -“3,; “5, = .o1 = -“5,,. 

aoe Stable? 
.2982x10 > no 
2163x107 no 
-1245x10°* no 
.2264x10 * no 

=:3403x10"" | __yes 








— 
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TABLE ie ab 








(confer) 














haxrgest 
Re 7 0) 


Peso 
.2090x10- 
.2018x10> 
.1950x10_ 


.1883x10 


.1996x10- 
.1595x10- 
1497510" 


.1499x10_ 


.2163x10- 
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i 


1 


al 


i 


p 


dh 


1 


A 


a 


ih 





Stable? 


no 


no 


no 


no 


no 


no 


no 


No 





TABLE IV-1 





(EGON 
k _k _ He - Pee) kK : _ =k 
i 1. =e 3y = 005 = 3 pi Sy, =a Ypi 
k nok k fae k ees 
an == Ani 4 a 4 i On = 6, 







Stable? 





Largest 
Rt. # 0 








or 


















soo esl’ no 
.1493x10°* no 
.1181x107+} no 
9429x107 no 
~- LOS 0scEG ie “no 
.6420x10° e no 
|) + 4783x10 | oe no 


Same, but il. = 





04 el 1449x10 MG yes 
04 se ~.2810x10 72 yes 


8 4 





>oO 
oO 
Oo 
© 
o 
QO 


Largest g 


Root x 104 





O> OnLy J), £0 
A=> Ae, Rig 0 





~2 40°! Tec pos, 5 


ky vs. Largest 


T Ome | At=0 





{120 

100 
iPareest of 
Root x 103 


60 


A0- 





40° 5 & 
Magnitude ka 


40° 3 5 ve 610 


K3 Wome? CSitmmnoOr 


Pome A To 





PUGURE IV~2 
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Waroes t 


3 eee 
Reot x"10- 


a0) 


aoee 





ay and ils VSuweencest Root 


.O1= 


Cissne- . O05=-k3,,5 k 


D/\ 





Ail At=0 ) 


FIGURE IV-3 


Le 
oe Orb, hi, VARY 
Largest on © : ADhy ds, VARY 
Root x 10° © 
18 
16 A 
Ls 
14 a 
OURS 4 2 3 4 5 a—h, SCALE 
a a gs p 
O .°4 .o2 o4 SOGIOT.0% 40 |), SCALE 






Ky and |kyl vs Largest Root 






(hy ,=Kyp=2.7 kg .=-.005=~kg, 


ko =eOl=-ks; All at=0 ) 
FIGURE IV~4 
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(iene emia -OoCrtt1Onal Constants for both ships, 


(Ky. k., and Keds set at values within their ranges, the rate 


3 


Bonstants, (k kK sea) Keds were systematically increased. 


ad 
This step is summarized in Table IV-1, part 5, and the effects 
rox Ky and Ky are illustrated in Figure IV-4. As can be seen, 


increasing k. decreased system instability, and increasing k 


ze 4 
had the same effect until Ik, | reached about .07, after which 
instability increased. Varying Ke had no effect on the 


largest root, though the other roots changed slowly as Ke 
increased. 

In part 6 of Table IV-1, the effect of increasing all 
three rate COllctales, Wlettehne proportional constants fixed, 
1s shown to be a decrease in instability. And in parts 7 and 
8, the same trend is shown with the Ky values being set at 
Mmiast 3.0, then 4.0. This last combination is seen to result 
in a stable system. 

Some investigations with non-zero time lags are summarized 
mn Table IV-2, and in Figure IV-5. As can be seen, in all 
cases investigated an increase in time lags caused the system 
to become increasingly unstable. 

ia Table IV-3, the results of giving only one ship sate 
Mmomatic Control system are summarized. For this case all 
io OU. Sthis Geable undicates that a stable system can be 
Echieved with only one ship automatically controlled, with 
control constants in the expected ranges. Physically, this 


Wn 


would appear to indicate that one ship has been "captured" 
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UCASE. LV-—2 


summary Of Stability Roots Obtained for Various 


Combinations of Control System Constants With Time Lags 


Cee Soe ye = 15 kts. 


Ship A LBP Bvemoe onto 8B LBP = 52 os. 


Puiliownda sted ks = 0 







k ~. af a k _ mae 
Ln x 1, = 2.; 27 = .01 = 2p 
At, At. At, Largest Stable; 
A A Re, =a 













2 .2399x10 “| no 

4 .2637x10 | no 

6 .3116x10 “| no 

Same, plus Ko 
f -U A 

-j 

22465510 no 

Oe econ no 

ool ae no 
~j 

ZO 7 lO no 

Poaees io) eno 

Fop7 x10 -) no 


8& 





feEOnm .) 


# <0 i> 0 Largest Stable? 
Leo teceleeRt. ~ 0 





. 2399x110, no 
P2637210 no 
Doe on no 

oe tens a k mks 

= -~,04 = eto = 1S Ge 

g 2 .1780x10 4 no 

g D .2362x10 | no 

8 Z .3473x10 7 no 








oo 








30 
Rmarcest ; 
© 
Root x 10° 
20 
O eae 
O 4 2 4 Rete le 
OG < 2 4 (6) <-At,, EOtrg 


Mame Wecs vs. Largest Root 


Cee ota ke = =, 00 Dace. & 


ke =,.0l1l==k 
ee 53 ? 





FIGURE 1V-5 


O10 





ee eo oe 


= am tet ei re cet ty 








ee ee 














ee 


Gamncryeot Stability Roots Obtained for Various 


Combinations of Single Ship Control With No Time Lags 


Ceo Pome — Ole Vv = 15 kte. 


oueee pe = 520.5 ; Ship B EBP = 528.5! 








SOAs Adie diet) 


eae eel 
























i: >0 Largest Stable? 
Rts. | Rt. #0 
2 aoe aon no 
2 3092810 - no 
2 .3914x10 > no 
0 Sa eels. yes 
eos oe Pies. es aieaea 
ea ee OC 0s, a Oe Se; 
i 
1p = 0 
k 
7 
ae 0 | 10 eee ae yes 
Ny 0 10 0 -.6003x10 > yes 
eee een eet Eee en See ea eee 
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Me ethe Other through the interaction force and moment. This 
was an unexpected result and it was felt that it might be 
indicative of an error in the formulation or solution of the 
equations. 

To investigate this possibility, some extra computer 
runs, for Special circumstances, were made. In these runs, 
ence MiccCmiciwen Berec, moicent and derivative inputs for 
program SOLVE were set equal to zero. This is equivalent to 
the ships being infinitely far apart, where no interaction 
Serectes Can be felt. in this circumstance, it should prove 
possible to obtain a stable system with reasonable control 
constant values for each ship, and, this, in fact, was the 
case. 

A further check was made, where the ships were given 
control systems with different constants---then run again 
with the same two sets of control constants, but switched 
from one ship to the other---both runs made with the inter- 
action effects all set at zero. These two runs should result 
in identical sets of stability roots if the solution is 
POcreecly @Mene, and, dad, in fact, do so. 

The last check run to verify the accuracy of the solution 
method was with all control constants zero for both ships. 
The interaction effects were all left at zero, as well. This, 
PaCrsect olen the Ships intinitelyetarvapart, with no 
automatic control, and should result in neutral directional 


mootiity,. ° Phe result of this run was four zero roots and 


22 





mix negative roots, indicating a dynamically stable systen, 


Meutrally stable, probably with regard to heading 


my. Conclusions 


The investigation reported herein has prompted the follow- 


mae conclusions: 


A. 


The methods of formulation and solution of the equations 
of motion appear to be correct. 

It appears possible to achieve a stable systen of 
iPiGmotips Ol parallel courses in close proximity, 

using automatic control. 

eu oMecuabMety io Ost art eCeted by comirol SenSitiva ty 


to heading angie error, ---that is, k, terms; and 


1 
least affected by sensitivity to longitudinal pos-~- 
titOmeernOoreand error rate, and ~~~-~that is, k « 
and Ke Tenens . 

it aoe possiple to achieve a stable system with pro- 
portional control only, (kp) ky, Ke Os 

Pinto mweesoiple, in isolated circumstances, at least, 
to achieve a stable system with only one ship 
Mpoving  alwromatic control. 

Under some circumstances, increasing the magnitude of 
a coutrol constant, and thereby increasing system 
twee cOeOne Of the €rror parameters, can have 


Seo cmeT oe ort eGt On SYStem Stability, (as illustrated 


iieercure LV¥-2), 


oS 





V. Recommendations 


The fuiliowing recommendations are made fur future 


investigation into the practicality of employing automatic 


Bontrols for ships engaged in operations alongside and close 


aboard, as in underway replenishment: 


A. 


A detailed examination of the theoretical development 
used in this thesis should be conducted by someone 
Ouner than the 2uchor, to remove any reservation as 
HOUT wacGliiacy ©: tne method - 

The examination of system stability should be made 

for two ships whose full set of hydrodynamic derivatives 
are: known. (+his would eliminate to need to 
Synthesize the eee ner ce of two different Siesps , 
as the author did with the OLNA/MARINER modification). 
iierrnesiiitsaor any Such theorectical inmgesti cation 
should be confirmed by model tests. 

The method should be extended to examine the effect of 
including en OULSide excitation of the system, Such 


eomviawics wild, tC. 
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Appendix 7. 


’ COMPUTER UTILIZATION 


jm  biscussion 

Three computer programs have been written as part of this 
investigation. The three programs, described in Appendices II, 
TII and IV, are used in sequence if the investigation of the 
Stability of a two shiv system, such as that studied herein, 
Starts at the level of known information with which this 
thesis was begun. If more information concerning the inter- 
action of the two ships being studied is available, it may be 
possible to enter the sequence at a later stage. (This 
Sequence of the programs is illustrated in the overall flow 
chart, Figure AI-1). Precisely because any possible future 
use of these programs may commence at a different stage, the 
three programs were not unified into a single large program. 
The use of the programs will be summarized in the reverse 
order from that in which they were used in this work. 

Program SOLVE, the third program used, is the heart of 

\ 

Pees COMputcer solution and actually solves for the system 
Stability roots. It also requires the greatest amount of 
Known information. If all the interaction and open water 
mydrodymamic derivatives for Ship ene the interaction 
derivatives for Ship A in Ship B's position, are known, 
program SOLVE may be used directly. 
¥ PeGeeen DETERM Can be used to gel the force and moment 


derivatives with respect to propeller speed, n, and the 


oo 





mreeraction derivatives x and Bhs for both ships, if these 
Em@entitbes abe net known. To do this, DETERM uses the values 
of Y and N at the positions of interest, plus You Ns and = 
Open water for Ship A as inputs. EHP and RPI data for Ship A 
are also required inputs. The output of DETERM then forms 
pet Of the required input for SOLVE. 

Program INTCOF is used to determine the interaction 
derivatives of Y and N with respect to a and £8, and the 
values of these quantities at the vositions being examined 
can then be used as part of the SOLVE input. In the vrocess 
of getting these inpecenion derivatives throughout the range 
of interest, program INTCOF also calculates Y and N everywhere 
in that range, uSing the interpolation technique described in 
Seapter Jil and illustrated in Figure III-l. 

The way in which INTCOF and DETERM were written was 
largely dictated by what information was available concerning 
the interaction phenomena and were steps toward getting all 
the derivatives needed for SOLVE. AS a result, INTCOF, 
especially, would hardly be likely to be useful to an 
investigator starting with a different level of known inter- 
Fecwvon information. lkteerOor Instance, Y and N were known 
from model tests everywhere in the range of interest, the 
interpolating portion of INT'COF would not be needed and the 
derivatives with respect to a and 8 could be obtained more 
Cieeetiy). Program SOLVE, however, should be directly useful 
for two ship system stability investigations for any two ships 


whose hydrodynamic derivatives are all known. 
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Appendix II 


PROGRAM INTCOF 


I. Discussion 
Program INTCOF takes the values of Y and N where known 
from reference A-1l, interpolates to get Y and N everywhere 


in the range of interest, then calculates val Yer Nu and N 


B 
everywhere in this range, using subroutines DERIV and DERIVB. 
Program INTCOF is highly specialized to the use of the 
data of reference A-1l. Any investigation using Y and N 
PreOorMatson from another source to get the above four inter-_ 
Pear TOnNacLivatives would be unable to use INTCOF. And, of 
course, an investigation using the [A-1] data could use the 
results herein, directly. For these reasons, the input format 
Bowereed £Or INTCOF will not be described in detail, and the 
remainder of this Appendix will merely contain a flow chart, 


Beogram listing and results for INTCOF. 
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I} Frow CHrarr —- PRocRam Iwrcor INTCOF PGL 
7 oe 
T 'rrom | 
\INTCOF | 
1 t r 
NRUN= 1 | ae Lay. VARIABLE 1S 
ba = 
mae Y(M,NJ= Y(MM N) + 


N-MOMENT OW 
PyNGew ASN) = 


SHIP AY 
Y(MM,N)) ® (CY (i, 
MIDALF)— Y(MM, 
MIDALF))/ Y(NBETAS, 


STAR 








NGETAS, NALFAS, 
(y(4,6), N= 4, VALFAS) 
(Y(NGETAS N), N= J, 














ee MIDALF )~ ¥(MM, 
MIDALF ))) Ke- 4, NALFAS 
MIDALF = . 
(NALFAS +1) /2 FROM INTCOF 
<0 Powe ts 
END a. . fae 
(Y (N, MIDALF), ee Do seo AL= ALINIT + (K- 1) ¥ 
N= 4, NBETAS ) 2 ALING 
? 
YES 
ALINIT, ALING, 
BEINIT, BEINC END 
DO 3800 P 
d 
NAENOs NALFAS -41 ? 6 


NBEND= NBETAS~4 YES | Je 2, NBETAS 


is 
Ma oy BEND 






a” 


BE = BEINIT + (J-1)¥ 
BEINC 


BE, yv(3,K) 


NRUN = 2 2 TRVE 





FALSE 









VARIABLE 3S 
Y FORCE ON 
Sty on 











/ a 
N<- 4, NALFAS 









NOY DO 


acre ec | —— 


FALSE 


TO INTCOF Yes 


Eee ss 


¢ 
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FROM INTCOF INTCOF PG, 2 
PG 28, #1.2 ae 


Tro IN Peor 


vy PG. 1, te 4.4 





yes 


NOBEND =NBETAS -2 


NDAEND = NALFAS -2 









= 
M< 4 NBETAS 
END 
NonoO - LOO 
> 
BE TA= BEINIT + ? 
(M-1) ¥ BEING YES 








ALFA=ALINIT + 
(M-1) *% ALINC 








DERIV (NALFAS, 
ALINIT, ALINC, M, 
Ver oLoce) 









NRUN=3 


#124 
emi 


DERIVB (NBETAS 
BEINIT, BEING, M 
Y, SLOPE) 


) 


— 1, NDAEND 


0 
ALF = ALINIT+L¥ALINC 


L<- 14, NOBEND 


EO FIN Te OF 


BETA, ALF, | ¢ ; BC aderat 415 
SIEOPE (tL) BET=GBEINI\T + 


ea BEING 








= 


OS 





DERIv PG. ji 
(INTCOF) 








Jusroutive Deriv (NXY, XINIT, XOELT,M,Y, SLOPE) 


FROM INTCOF 
PG. 2, # 2.3 







K = 4 
X= Yel T 


REP = + 2 XDELT 
KP=K+2 
SLOPE (K)= 0.0 


SLOPE (K) =(¥(M_ KP) 
—Y¥(MK))7 
UP = xX) 





AMID = XAADELT 
K = oe 
Y= KMID 


KP =NXY i FALSE 
TRUE 


RETURN 
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DERIVE PG. 4 
ther Coe) 
SuBroutine Derive (NPT S BINIT, BINSM,Y, SLOPE ) 
FROM INTCOF 
er cme 1c.4 


V 


K= 
B= BINIT 
















BUP = B+ 2.¥ BINC 
KP=K+2 
SLOPE (K)= 0,0 
‘SLOPE (K)= (Y(KP,m) 
-Y(K.M))/(BUP -B) 


BMID= B-+-BINC 
K= Kt 
B= BMID 





V 


RETURN 
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Appendix III 


PROGRAM DETERIM 


aie Discussion 

i eee oe, OUtpuc Format 
eee wronvecnart 

i esaeeawan Listing 


V. Sample Results 





mependix £1 Lf 
PROGRAM DETERM 


ieee Discussion 

Program DETERM determines the hydrodynamic derivatives 
with respect to propeller speed, n, which were not available 
in the literature surveyed. It also determines the variation 
Bemeoeecaction Cttfects with velocity, u. Thus the outputs 
arc X, Yn and No Ore ona A; yy and Ny for Sime, A; ame Ya 
and Ny Lenser ia oie 2 S DOSition. These latter values, 
when used aS inputs in program SOLVE are modified to their 
Pppeopriate values for Ship B in Ship B's position. (xo Yn 
and No are position-independent. ) 

The values of Ree and No are calculated by using the 
known values of X, Y and N exerted by the propeller at - 
Seradldbedum speed in open water and assuming that propeller 
forces ana mements vary as the square of propeller spced. 

The derivatives with respect to u also are calculated by 
taking the known values of Y and N at equilibrium and assuming 
that they vary as the square of the velocity, u. (The 


procedure followed in calculating these quantities is 


described in detail in Chapter ITI, PORCEDURES. ) 


Loo 





ee Input/Ourput Format 


Program DETERM's input is ordered as follows: 


LITEM 


LBP Ship A 


Beam Ship A 


peewee at 1 kt. 
equil. 


reel eat 
Eee at 
au 1 . 
Per at 


EHP at 


Beeed at | kt. 
equil. 


~~ 


RPM at 


EHP at this speed 


# positions examined 


Long. 


Lat. separation at 
each position 


Y for@e at each 


this speed 
this speed 


speed (kts.) 


equil. 


equil. 


this speed 


separation at 
each position 


WoSstils On , 


MN moment at each 


pPesition, 


specd 


syeed 


PROGRAM 
SYMBOL __ 


7a P 


ABEAM 


weer) 
RPM (1) 
BHP (1) 
ae) 
RPM (2) 


iAP (2) 


V(3) 
RPM (3) 
EHP (3) 


NPOSIT 


ALFAO (TI) 


Soon a ) 


CARD 


6 


COLUMNS 


SS me 


a= 20 


€ard for each position) 


AY (iM) 


AN (M) 


eS, 


7 


7 


J, -2L4 


T5s20 


(one card for each position) 


FORMAT 





PaO. 2 


Halo. 2 


PLO. 2 
Or 
BeeO re 
PEO rez. 
PO 


FLO. Z 


F10.2 
F10.2 
Ue 


T4 
AO 


ELOe 2 


F1l4.4 


F1l4.4 





LTEM PROGRAM CARD COLUMNS FORMAT 


ee 
merorce on Ship A 
in Sip B poOSLtion BY (N) g i 14 E14.4 
N moment on Ship A 
in Ship B position BN (N) 8 oe ba. 4 


(one card for each position) 


Ys Ship A AYDEIL 9 es B14.4 
Ns Ship A | ANDEL g oS 2 B14.4 
Cmaeequil, Ship A ADELO o oe oe One 


The output format for DETERM is illustrated by the 
Sample results section of this Appendix and should be self- 
explanatory. In the output, XO is represented by DX/DN, 


A i DY/DN, etc. 
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— das 
oe 

NPOSIT, 

fe. FAOWT) STSO(T} 

et, NPOSIT) 

(AY (mM), AN UM), 

M= 1, NPoSIT) 
(BY(N), BNON} 
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ADEL © 
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40 
ee - 1S ) 
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Ca 


X= EMP(M)¥ 550. / 
(V (4) ¥ 16878) 


ee wee 











PROGR AM Dererm 


o. ics 


AXN= (( X13) ~ X(4)) / 
(RPS(3)— RPS(AI))/ 
(0.5 ¥RHO XHALRPXX3.) 

% V(2) * 1.6878) 


‘2 


DAYD= AYD * 0.5 XRHO 
¥ (ALB PXK 2.) ¥ (VO2) 
RAGE (BIH 2. 
XADELO/57.3 











DAYDUP= DAYD x 
(RES (3)/ RES(2)) ¥¥2. 





















SLAYD= (DAYOUP — 
DAYD) /(RPS(3)—- RPS) 
AYN= ~ SLAND/ (0.5% 
RHO X(ALG PX % 3.) ¥ 

VIZ) LGB7B) 


DAND= AND %0.5 * 
RHO x(ALS PRS.) x 
((V(2) 49.6976 )xX*2.) 
¥ ADELO/ 7,3 
DANDUP= DAND * 
(RES(3)/ RPS (2) #82. 
SLAND = (DANDUP~DAND 
ARPS(3)- RES (27) 
ANN=— SLAND /(0.5 
YRIO X(ALBPH%4,) 
¥V(2)% 46678) 


FROW 
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DETERM 
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¥XZ.)XALOPH*X2.) 
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Vi) ¥¥2, 
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AYU= SLAY/(0.5#RiO 
#&V(2)*X1.6878 * 
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PG, 2, #21 
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cal tt 4.4 
DAY= AY(L) ¥(0.5 ¥ 





FROW DETERM 


PG. 9 et 
7 


DAN= AN (L)X (0.5% 
RHo x ((v(2) *4.687 8) 
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DANUP = DAN (v(3)/ 
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STSD(L), 
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‘FOR SHIP A 
IN SHIP B 
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DBY= BY(L)* (0,5 RHO 
%((V(2)%4.6978)%% 2.) 
*ALB PY 2,) 
DBYUP= DBYx (V(3)/ 
V(2)) are. 
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RHO x ((V(2)%4,6878) 
VED PKXALBP XY 3S, ) 
DSNUP = DBN x (vis) / 
V(Z))#8 2. 

SLON= (DRNUP-DBN)/ 
((VI3)-V(2)) ¥1.6876) 
BNUINT = SLBN/(0.5 
RHOXVIZ) X41, 6078 
(ALBPx** 3.)) 
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Appendix IV 
PROGRAM SOLVE 


im Discussion 

Program SOLVE, using the hydrodynamic derivative inputs 
for Ship A, and for Ship A in Ship B's position, calculates 
meemtnll set of Gerivatives for Ship B. The derivatives for 
poh ae are dimensionalized in the program, and the 
applicable coefficients for the variables of the equations 
of motion are computed. SOLVE then calls subroutine DETER 
yhich computes the values of the coefficients of the polynomial 
in) which is the result of evaluating the determinant of 
eeenwequations of motion coefficients. 

Subroutine DETER, in the process of performing necessary 
molynomial addition, subtraction and multiplication, calls on 
[eerary yOutines of the IBM Scientific Subroutine Package (SSP). 
These subroutines are PADD, PSUB and PMPY and are fully 
Beseribed in the IBM Apolication Program publication H20-0205-3, 
entitled "System/360 Scientific Subroutine Package, (360A-CM-03X) 
Version JII". These routines are included in the program 
listings in this appendix and were used without modification. 

Once DETER has evaluated the determinant to get the 
Gwemacteristic polynomial in oe execution is returned to the 
ime program, SOGMVE, sto determine the roots Of this equation, 
which are the fetertlity Boots Wor the system. ~~ ine roots are 
determined by a fourth SSP subroutine, POLRT, used without 


modification and listed in this appendix. 
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ee ee meme TS Ae rk gm Sl cl 


mi. inpur/Oultout Format 


Program SOLVE's input is ordered as follows: 


LTEM PROGRAM CARD COLUMNS FORMAT 
3 SYMBOL _ 
phip A Characteristics ALBP, ABEAM, 1 70 7PLOew 


ADRAFT, ADISPL, 
ACB, ACP, ACM 


Ship B Length BLBP 2 ke PLOeZ 
Parl. Spd. in kts. VO 3 eee) LO 
Dimensionless open 2 AXUD, “AYU, eae4 B14.4 
water derivatives, awe, AYPD, AYPOD, 
Ship A Peet, ANV, ANVD, 
(AXU="u, , ee) ANPD, ANPOD, ANDEIL, 
A= 1-5 
Dimensionless AO AYN, 16-18 1-14 AeA 
derivatives from ANN 
Pie RM 
# Positions examined NPOSIT io 1-4 T4 
ffong'’l and Lat. ALFA(I), 20 i220 2F10.2 
Separation at each sD (I ) 


pos 1tion 


(one card for each position) 


Dimensionless vol Peet), AYB{E), ee an ae 4 
Y,, N., N,, Ship A ANA(I), ANB(T) 
B O 8 21 


{one card for each position) 


bimensionless Y , Ay (iL), 22 L263 2 aa 
es ANU (I) | 
Nw Ship A 


(one card for each position) 


> 


Dimensionless Ye Bo tie (1), 23 1-28 2h ee 
Nae Ship A, shia: eat) 
Position 


(one card for each position) 
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ITEM PROGRAM CARD COLUMNS FORMAT 


_ SYMBOL on 
Pamensionless Yue Yor BYAUNE (I), 24 Le 2E14.4 
Ship A, Ship B eee (1) 
position 
(one card for each position) 
Dimensionless Noe Nee BNAINT (I), Do le 28 2E14.4 


: 
pip A, Ship B a St) 
position 


(one card for each position) 


The output format for SOLVE is illustrated by the sample 
results section of this Appendix. In the output, vay 1s 


represented by DY/DA, etc. 





Te. SOLVE Flow cHARY 


Procram Sotve 
START = 


| % 
AYA(I), AYB(T) 
ANA CT) ANG(T), 
AYV(T), ANU IL) 
BYUINT (I), BNUINT (1) 
BYAINT(Z)) BYGINT (T) 
BNAINT (I) BNBINT (T) 
owe, NPS IT 



























ALBP, ABEAM 
ADRAFT, ADISPL, ACB, 
ACP, ACM, BLBP 













BBEAM= ABEAM x 
(gLbe /ALSP) 
BRRAF T= ADRAFT * 
(pipe /ALBP) 
BDISPL= ACB ¥ BLBPx 
BBEAM *BDRAFT/35. 










DIM1= 0.5 ¥ RHO * 
ALBP Xx 3 





ALSP, BLBP, ABEAM, 
BBEAM, ADRAFT, 
JBDRAFT, ADIS PL, 
BDISPL, ACB, ACB Ace 
ACR ACM, ACM DIMG = DIMS ¥ ALBP 
DIM7=O0.5¥ RHO | 
¥ALBP XUOKX2 
DIMS = DIM? KALB P 


|UO= VO ¥ FPS DIMS = DIMS *ALBP 


AXUD= AXUD *DIM4 
AYVD = AYVD *DIM4 
AYPOD= AYPDD ¥DIM2 
ANVD = ANVD ¥ DIM2 
ANPDN=ANPDD + DIM3 
AXU = AXU ¥ DIM4 
AYV =AYV ¥ OIM4 
AYPD= AYPD ¥DIMS 
ANV= ANV ®DIM5 
AXN= AXN ¥ OIMS 
ANN = AYN*¥ DIMS 
ANPp= ANPD ¥ DIMS 
ANN = ANN *® DIM6 


a. 


AXU, AXUD. AYV, 
AYVD, AYPD, AYPDD, 
AYDEL, ANV, ANVD 











ALFA(L) STSDIT), 
I=1, NPOSIT 





Cae 


1590 


a. a, 
V 
AYDEL= AYOEL x DINM® 
ANDEL= ANDEL *DIMS 
Ci= BLBP/ALGP 
Cie: Cy} KA? 
C3=C1¥XZ 
=—_ C4 C2 x*z2 
ee CS = Ciges 5 
yUr Axum 
DIMZ= DIML ¥ ALG P b oe 
BXUD= AXUDXC3 
DIM3 = OIM2 ¥ALBP 
BXN-= AXWN¥ C3 
DIM4= 0.5 ¥RHO %XU0% BYN = AYN * C3 
ALBPx*¥2. BYDEL= AYDEL XC. 
DIMS = DIM4* ALBP Byv= AYV¥C2 


SOIVE FG4 







BYVD= AYVD %C3 
BYPD= AYPD * C3 
BYPDD= AYPDD *#C4 
BNN= ANN Cd 
BNDEL.= ANDEL *¥C3 
BNV= ANV ¥CS 
BNVD= ANVD *¥ C4 
BNPD= ANPD xc¢ 
BNPDD= ANPDD *CS 













AXY, BXU, AXVD, BUD, AV 
BYV, AyvD, GyvD, AyPD, 
BYPD AYPDO, BYPDD, 
AYDEL, BYDEL, ANV, BNV, 
ANVD, BNVD, ANPO,BNPD 
ANPD)D BNPDO, ANDEL, 
BNDEL, AXN, BXN, ANN, 
BYN, ANN, 8NN 


TO SOLVE 
PG. 2, #24 





FROM SOLVE 
Ce 


’ 
( Ant 
H2.4 


1000 
oe deers | 7 











2tZ.4 
AYACK)= AYALK)4 DIMT 


AYG(KY= AYE (K)¥ DIM 7 
AYU (K)= AYU (KK) ® DIM 4 
ANU(K)= ANU(K) *OIMS 
AWA (K\= ANA (K) ¥ DIMS 
ANB (K)= ANG (K) XOINS 













BYA(K) = BYAIN ST (K)? 
C1i*DIM7 
BYB(K)= BYBINT (K) # 
Ci* DIM? 
BWA (K)= BNAINT (K) ¥ 
C2*DIM& 
BNS(KY= BNBINT (K) % 
C2 ¥ DIM 8 
BYU (icy = BYUINT (K)* 
C2% DIMA 
BNU(K)= BNUINT (K) ¥ 
C3 * UiM5 





— —- 


BETA(K)= STSD(K) + 
(ABEAM+ BREAM) /2 


| 
ALFAIK) STSDIM, 


BETA (K) 





















AYA(K), BYALK), 
AYB (K) BYBIK), 
ANA(K), BNA(K) 
ANB (K), BNB(K), 
AYU(K) BY UTK) 
ANU(K) BNU(K) 


1004 
MI<— 1, NCONS 
















ARTA, ADT2, 
pod, BDT2 












(AKIM) M= 14,6), 
(GK(M) M= 1,6), 
ADT1, ADT 2, BDT4, 
BDT2 


ilo 






SOLVE OG 2 




















Axit4)= AXW ® As) 
j 


ETC. 
DETERMINING POS{|TION- 


INDEPENDENT, CONTROL- 
DEPENDENT CoEFFIC-— 
(ENTS eCr pO re 


VARIABLES 


f 
i 
! 


QNG(3)=BNPOD — BNDEL 
BK(2) ¥ BOTY 


ie 


AYZ (1) = AYB(K) + AYDEL* 
AK(3) 
AN2(1)=ANS(K)-+ ANDEL ¥ 
AK (3) 
BY2Z(1)= BYB(K)+ BYOEL * 
BK (3) 
BN2(1)= BNS(KI+ BNDEL> 















DETER (ToET) 






_—_4t2.3 


ora me eke 


TO SOLVE 
PG 3, een 





Prom SOLVE 
P6.2,  #2.3 









oe 2. 0. 
Sooo 


ame LDET 


NOS = et 11) / 


POLRT ( ADET, COF, 
TORD ROOTR, 
ROOTL, TER 


ae 


TER, (ROOTRIM), 
M=1,IORD) 
TO SOLVE 


Eo PG, ic /\ ae 2. 


e522 
END 


DO Vege 
160+ 






TO SoLYE 
OG 2 #24 





SOLVE 


PG. 
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